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1.0 SUMMARY 


This document describes the computer program L216 (DUBFLX), and contains the 
following: 

o Summary of analytical development 

© Outline of the program structure 

© Specification of information necessary for execution of the program 
o Control cards 
0 Required resources 
o Input data 

® Description of program output 

® Listing of program restrictions 

Q Listing of error diagnostics 

o Description of sample problem 

Program L216 computes three-dimensional, subsonic, unsteady aerodynamic 
characteristics for arbitrary configurations that can be modeled as combinations of 
lifting surfaces and slender bodies. Input to the program consists of configuration 
geometry and aerodynamic condition specifications read from cards, and modal data, 
which may be either read directly from cards or interpolated from arrays of modal 
coefficients read from an externally developed magnetic file. 

Computations are based on a finite constant pressure panel concept, the doublet lattice 
method. The primary function of the program is to evaluate the kernel expression in the 
integral equation relating pressure and normalwash on lifting surfaces. If modal data is 
supplied, the geometric normalwash may be determined and pressures and generalized 
air forces calculated. The pressures may be integrated with appropriate weighting 
factors to produce normal forces, moments, and stability derivatives. 

Optionally the user may file for subsequent utilization, geometric data, matrices of 
normalwash factors (used in aerodynamic influence coefficient-type calculations), and/or 
pressures and generalized force arrays. 



2.0 INTRODUCTION 


Computer program L216 (DUBFLX) may be used as either a standalone program or a 
module of the program system DYLOFLEX (ref. 1). Program L216 is a modification, 
accomplished under contract NASI -13918 to satisfy contract requirements developed m 
reference 2, of the program described in references. The tVieoreT,ifT.t 1 '.trniu t ut j< )H ntid 
tlie eq- ration solution Tiietho>l is that of refevesce 


Modifications necessary to satisfy the requirements of DYLOFLEX included: 
« Addition of modal interpolation routines 
9 Modification of input and output (contents and formats) 

• Improvement of internal documentation (commenting) 
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3.0 SYMBOLS AND NOMENCLATURE 

The following list contains items that appear in this document except section 6.3 (card 
input). 

Engineering 

notation Definition 

2 

A Reference area (length ). 


c,Cj. Local, reference chord length. 

Cn, Cm Section normal force, moment coefficients. 


Cz, Gy 

^mi 

ACp 

ACz, ACy 
^rs 


Normal, aide force coefficients. 

Pitching, yawing, rolling moment coefficients. 

Pressure difference coefficient. 

Slender body pressure coefficients. 

Element of matrix of normalwash factors-lifting surface on lifting 
surface. 


e Element semi width. 

F|.g Element of matrix of normalwash factors-slender body on lifting 

surface. 


f 


Mode shape. 


f 

g 


Slender body off the plane-of-symmetry factor. 


h 


Modal displacement. 


VT. 


k Reduced frequency (o>Cr/VT ) (radians). 

K(x,y,x,o>,M) Kernel function. 

M Mach number, 

q Generalized coordinate. 

Q Generalized force. 

Rq(x) Radius of slender body (length). 
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s 

Vt 

w(x,y,z) 

Wg(x,y,z) 

x,y,z 

Ax 

y 

6 


a 


T 


4 >' 


Subscripts 

a 

B 

i 

j,k,8 

le 

r,s 

R 


Reference semispan (length). 

Aircraft velocity (length/sec). 

Normalwash. 

Normalwash induced by slender bodies. 

Augmented normalwash. 

Surface coordinates; receiving point coordinates. 

Element length. 

Element dihedral angle (radian). 

Symmetry factor (Y = 0). 

Symmetry factor (Z = 0). 

Element coordinates; sending point coordinates. 

Streamwise coordinate. 

Spanwise coordinate. 

Radial coordinate. 

Modal translational deflection. 

Streamwise derivative of the modal translational deflection 

d4>/dx. 

Axial line element. 

Slender body. 

Modal degree of freedom. 

Strip, box, line doublet elements. 

Leading edge. 

Receiving, sending elements. 

Reference. 



Superscripts 


(f) 


Slender body. 


Differentiation, d/dx. 

Matrix Symbols 


u 

[] 

n 


Diagonal ma.trix. 
Rectangular matrix. 


Column matrix. 



4.0 ENGINEERING DESCRIPTION 


A complete theoretical development is contained in reference 3 Pertinent equations and 
solution methods are discussed in this document in section ^ • 1 and applications in 
section 4.2. 

4.1 SUMMARY OF THEORETICAL DEVELOPMENT 
4.1.1 EQUATIONS 

Consider a linearized representation of three-dimensional subsonic unsteady 
compressible flow. Given these assumptions, the velocity normal to an oscillating 
surface is related to the lifting pressure by the integral equation: 


w(x,y,z) =-^ // K(x-^,y-i 7 ,z-^,a>,M) ACp d^ da 

OTT 

LS 


( 1 ) 


Discretizing the expression (i.e., assuming that the lifting surfaces can be approximated 
by segments of planes that are divided into small trapezoidal elements), the integral 
equation may be approximated as: 


w(x,y,z) = -^SACp // K(x-^,y-Tj,z-C,&j,M) d^ da 

OTT a S _ 


( 2 ) 


A further approximation in the doublet lattice method is that the integration of the 
kernel expression in the streamwise direction is accomplished by lumping the effect at 
the 1/4 chord of each element: 


w(x,y,z) =-^ X ACp J K(X -^1 ,y-T 7 ,z-S,eo,M) da 
OTT a ^ ® /4 


(3) 


A set of linear algebraic equations may be formed from the relationships above and 
written in matrix form: 


{w} = [D] {ACp} 


where a typical element of [D] is: 


Drs =4^/ K(x-^i/ 4 ,y-T 7 ,z-C,o>,M) da 


Stt 


(4) 
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If slender body interference effects are to be considered, the incremental normalwash on 
lifting surfaces in the presence of slender bodies is calculated using Woodward’s method 
of interference surfaces to determine lifting surface/slender body interaction effects. 
Proceeding analogously as with the preceding equations: 

WB(x,y,z) = ^2 ACzgRo8A^sKZg(x-fag.y-T1a.z-Sa.<".M) 


(5) 




AC Y^ROg A^g KVg (x-^ag .y-'Ha ,z-Ca 


This system of equations may be written in matrix form: 

{wb} = [Fz]{ACz} + [Fy]{ACy} 
where typical elements of the [F] are: 

Frs ~ ^I^Og^^s^rs 

where ACy and ACz are components ofACp^^\ determined using slender body theory: 

ACp^^^ - 2ir(R'ow^^^ + w^^^'Rq/Z 4- ikf w^^^Ro/c,.) (7) 

Then, for lifting surfaces with the effect of axial singularities introduced: 

{wj^} = {w-wb} = [D]{ACp} (8) 

Equation (8) is the form used to determine pressure differences on lifting surfaces 
considering the interference effects of slender bodies. 

4.1.2 SOLUTION METHOD 

Procedure 

1. Matrices of normalwash factors [D] and [F] are determined using equations (4) 
and (6), respectively. 

2. Geometric normalwash for slender bodies is specified and ACp^^^ determined using 
equation (7). 

3. Normalwash induced on lifting surfaces by slender bodies is calculated using 
equation (5). 
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4. Geometric normalwash for lifting surfaces is specified and combined with induced 
normalwash to form the augmented normalwash on lifting surfaces, Wjj. 

5. Pressures on lifting surfaces are obtained by direct solution of equation (8). 

Evaluation of the Kernel Function Expression 

Landahl’s expression for a nonplanar kernel function is used. In order to achieve 
increased accuracy, the kernel is separated into two parts: the steady component which 
may be integrated directly, and the unsteady increment which is evaluated 
approximately by fitting the kernel expression with a parabolic function and 
integrating analytically. In steady flow, the line of doublets is equivalent to a horseshoe 
vortex whose bound portion lies on the 1/4-chord line of the element. 


Normalwash Boundary Condition 


The boundary condition is merely the substantial derivative of the surface deformation, 
and for harmonic motion: 

{w} = - 

where 




{q} 


(9) 


[<^] = mode shapes which define translational displacements in the y or z 
direction 

[(^'] = the streamwise derivative of the mode shapes = d/dx [(/)] 


The displacement sign convention used in DUBFLX is shown in figure 1 





Figure 1 .—Displacement Sign Convention for DUBFLX 
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For liftings surface elements, the translational displacement used in equation (9) must 
be normal to the surface. jSlender body translational displacements can be either in the 
z or y direction depending upon the type of doublets associated with the slender body. 
The normalwash control points for slender body line elements are the midpoint 
locations; numerical experimentation indicates that the Kutta condition will be satisfied 
for lifting surfaces if the normalwash control point is chosen as the element midspan 
3/4-chord point. 

A slender body does not induce normalwash on the particular interference lifting 
surface associated with itself. An interference lifting surface does not produce lift due to 
its own motion; only the normalwash induced on these surfaces by other aerodynamic 
elements produces lift. 

When using doublet lattice in the DYLOFLEX system, the modal data used in 
equation (9) may be input in either one of three ways: 

• Defining a set of polynomials 

• Using modal interpolation arrays 

• Tabulating the displacements and slopes on cards 

The polynomials define the displacements as functions of the x and y coordinate of the 
surface. They are described by the user inputting on cards the order of the polynomial 
and the coefficients of each term of the polynomial. The modal interpolation (SA) arrays 
are generated by L215(INTERP) (ref.U )• These arrays are read from magnetic file and 
used with the OUBFLX geometry data to determine the needed displacements and 
slopes. Finally, the tabular modal input consists of inputting via cards the 
displacements and slopes at each control point. In all cases, the modal displacements 
must be defined consistent with the sign convention shown in figure 1. 

When using DUBFLX independent from the DYLOFLEX system, an alternate method 
does exist in L216 to describe the mode shapes. The doublet lattice program has a set of 
internally defined modes which the user may choose to use. However, the modes are 
inconsistent with the DYLOFLEX sign convention and are only defined for box control 
points. 


4.2 APPLICATION 
4.2.1 MODELING IDEALIZATIONS 

Arbitrary configurations are idealized as primary lifting surfaces (wing, tail, nacelle 
duct, control surface, etc.), slender bodies (fuselage, store, etc.), and, associated with 
individual slender bodies, interference lifting surfaces (see fig. 2). 

Primary lifting surfaces are represented by a series of planar segments, each of which is 
a collection of small trapezoidal elements arranged in strips parallel to the free stream; 
the pressure is constant over each element, and the resultant force is assumed to act at 
the midspan 1/4-chord location. Slender bodies, axially symmetric bodies with linearly 
varying radii, are represented with a system of line doublets located on the centerline of 
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Figure 2 _ Examples of Aerodynamic Modeling 
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revolution; the resultant force is assumed to act at the midpoint location. Interference 
lifting surfaces are considered as surfaces of revolution with constant cross-sectional 
shapes discretized just as primary lifting surfaces are. 

A comprehensive list of rules governing the discretization of the basic aerodynamic 
components (lifting surfaces, slender bodies, and interference surfaces) is presented on 
pages 49 and 50 of reference 3 • A brief summary follows. 

1. Lifting surface trapezoidal elements (boxes) are arranged in strips paralled to the 
freestream. 

2. Aspect ratio of the boxes should not be large; for the unsteady case an aspect ratio 
of order unity is recommended. 

3. Surface intersections and edges, fold lines, and hinge lines should coincide with 
box boundaries. 

4. Streamwise box length should be small relative to the basic wave length (see 
table 1). 

5. Boxes should be concentrated in regions where span loading changes rapidly or 
where normalwash boundary conditions are discontinuous. 

6. For planar or almost-planar wing-tail problems, spanwise boundaries on the tail 
must align with those on the wing. 

7. The radial distance between a normalwash control point and a box edge may not 
approach zero. 


Table 1.— Chord Division— Reduced Frequency Relationship 


Reduced frequency 

Chord divisions 

0.5 

4 

1.0 

8 


Note: Chord divisions based on relationship 



4.2.2 LIMITATIONS 


The doublet lattice method is subject to the limitations (common to other subsonic 
lifting surface methods) involved in formulation of the basic theoretical aerodynamic 
equations and of the kernel function expression. Accordingly, with these limitations, the 
theoretical modeling is restricted to a- linear representation of inviscid subsonic flow, 
and behavior of the kernel function deteriorates at high subsonic Mach numbers as flow 
approaches the transonic range. 

In addition, within this context two further approximations are made in the doublet 
lattice method. One is concerned with the approximate evaluation of the kernel function 
using a parabolic curve fit; this procedure, utilizing only a second order function, is less 
exact than that developed in the kernel function methods using higher order functions. 
It appears that this is not a major limitation. More significant is the assumption that an 
aerodynamic surface can be represented by a system of line segments of acceleration 
potential doublets for the purpose of calculating lift distributions. Practical problems 
associated with application of the doublet lattice method primarily involve this 
approximation; i.e., it is apparent that aircraft idealization is a critical factor governing 
the validity and accuracy of computed results. 
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5.0 PROGRAM STRUCTURE AND DESCRIPTION 


Program L216 (DUBFLX) is constructed as an overlay system consisting of a main 
overlay and eight primary overlays. 

(L216,0,0) L216 
(L216,l,0) INITIAL 
(L216,2,0) PRT2 
(L216,3,0) GEWSSL 
(L216,4,0) GEQ 
(L216,5,0) AMSOL 
(L216,6,0) AICSUB 
(L216,7,0) GIN 
(L216,8,0) SAMODE 

The main overlay, program L216, performs file addition and deletion, and calls the 
primary overlays as requested. Primary overlays INITIAL, PRT2, GEQ, AMSOL, GIN, 
and SAMODE are called for the typical problem. Primary overlays INITIAL, PRT2, 
GEWSSL, AMSOL, and AICSUB are called for the AIC problem. 

This program requires subroutines that are not part of the L216 code. Some are 
automatically obtained from the standard FORTRAN library when the- program is 
loaded. Others are stored in the DYLOFLEX alternate subroutine library, which must 
be declared at the time of loading. When the program is run as a standalone program, it 
is the responsibility of the user to generate input data in the format required by L216. 
Schematics of the L216 overlay structure and file communcations are presented in 
figures 3 and 4. The following is brief description of each overlay. 



Figure 3.— Overlay Structure of Program L216 
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NGETP NAETP 


Figure 4.— External Data File Communication 


5.1 PROGRAM L216 

This program intializes and deletes scratch files, reads input data, and calls the 
required primary overlays for solution and output data generation. The following 
keywords are used; the underlined characters are required. 

$DUBLA T Indicates that input data following is for use by L216 
$TITLE Printed and ignored 

$QUIT Terminates execution 
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S.2 PROGRAM INITIAL 


This program reads card input data for geometric and modal definitions, computes 
detailed geometric data, processes modal input, and stores geometric data for 
subsequent use in DYLOFLEX. The following keywords are used. 


REDUC ED FREQUENCIES 

GEOME TRY 

PANEL 

BODY 

MODES 


Introduces k*value array 
Introduces geometric data set 
Defines lifting surface geometric/modal data 
DeHnes slender body geometric/modal data 
Indicates modal data set 


5.3 PROGRAM PRT2 


This program evaluates the kernel function to form the matrix of normalwash factors, 
[D], for lifting surface elements, modifies the factors for symmetry, and writes the 
matrix of normalwash factors on scratch tape for later use. The kernel function 
expression is defined in reference 4 (appendices A, B and C). 

5.4 PROGRAM GEWSSL 

This program computes normalwash column vectors, { for each of the 

preselected AlC modes, combines these column vectors with the matrix of normalwash 
factors, [D], to form the augmented matrix of normalwash factors, [D| and 

writes the matrix on scratch tape for subsequent use. The preselected AIC modes option 
is described in reference 4 (sections 2.5, 2.6, and appendix D). The program generates 
AIC’s for lifting surfaces only; slender body elements are ignored. 

5.5 PROGRAM GEQ 

This program computes the normalwash matrix for all non-AlC analyses. If slender 
bodies are considered, matrix of induced normalwash factors, [F] , slender body pressure 
coefBcients, ACp , and the induced normalwash, Awq, are calculated. The augmented 
normalwash, [wfi], is determined and the augmented matrix of normalwash factors 
[D|wr] formed. The augmented normalwash and normalwash factors matrices are 
stored for subsequent use. 


5.6 PROGRAM AMSOL 

This program uses the quasi-inverse method to solve the simultaneous set of equations 
for the pressure on lifting surface elements. 

5.7 PROGRAM AICSUB 

This program calculates and prints the stability derivatives and pressure coefficients for 
the AIC option. Additionally, for this option, the program computes and prints harmonic 
gust coefficients. Analytic development of the gust matrix is described in reference 4, 
section 2.5. 
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5.8 PROGRAM GIN 

This program computes and prints the pressure coefficients, section normal force and 
moment coefficients, total lift and moment coefficients, and generalized forces. 
Aerodynamic results required for later use in the DYLOFLEX system are stored on file. 

Section coefficients on strip j per mode i: 


Cn.. = — .Axk 

!C; k 


-m 




( 10 ) 


Force coefficients on body f per mode i: 

Cy® =-||Ro,Ax,ACpj Y 
Total force and moment coefficients per mode i: 


Cz = (1 + S) [-W 2ej cos yjCjCn-. + 2 Cz/^^] 

i p A 


( 11 ) 


Cy. = (1 - 8) TCjCn-. + 2 g^^^ Cy/^^] 

1 A . jl ^ 1 


■1 


V r 2 / O 

Cm; = 1 Lc cmjj - ccnjj (xiej - xr)J 2ej cos yj 

Acr j 

-2 g^^^ I 


n. = - ccnji (xiej * XR>] 2ej sin yj 

+ 2 g^^^l2RopAxgACp ^ (xie- - xr) 

(f) g K Ki J 

Cg. = - 2 cCnji 2ej (yj cos yj + zj sin yj) 

a. V (r ^ r , 7 

+ 2 g (Cz; yc + Cy Ac ) 

(0 


(12) 
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Generalized force coefficients (half airplane): 


Qik 

S k 


+ 2g‘"S2ACp ®fjiAx»Ro. 
f J * * * 

where: 

(f) 

g = 1 -Bodies on the plane of symmetry 

= 2-Bodies off the plane of symmetry 


(13) 


The generalized aerodynamic forces defined by equation (13) are calculated on the right 
hand side of the equations of motion as defined in L217 (EOM) (ref. 5). 


5.9 PROGRAM SAMODE 


This program uses the SA (interpolation) arrays generated by L215 (INTRFP) to 
determine required modal data. For each lifting surface and slender body, the program 
reads the SA array, interpolates, calculates the modal displacements and slopes, and 
computes the generalized force integration matrix [BQ]. The modal data and 
integration matrix are stored on a random access file for subsequent use. 

The following is the element of integration matrix (box k/line elements, mode i). 


B^i — ^g(2e^ Ax^ 




Ax 




kki — s<f>kiqi 


s = 1.0 


(14) 


17 


6.0 COMPUTER PROGRAM USAGE 


The program was designed for use on the CDC 6600. The machine requirements to 
execute L216 are: 


Card reader 
Printer 

Disk storage 


Tape drive 


To read control cards and card input data 

To print standard output information, optional intermediate 
calculations, and diagnostic messages 

All sequential magnetic files not specifically defined as magnetic 
tapes are assumed to be disk files. The scratch random access file 
is always a disk file. 

For "permanent” storage of data, magnetic files are copied to and 
from magnetic tapes with control cards before and after program 
execution 


The program L216 is written in FORTRAN IV and complied with the "FTN” compiler. 
L216 may be executed on the KRONOS 2.1 operating system. 


6.1 CONTROL CARDS 

The following list is a typical set of control cards used to execute L216 using the 
absolute binaries from the program’s master tape. 

Job card 
Account card 




REQUEST(MASTER, F=I,LB=KL,VSN=6 6XXXX) 
REWIND(MASTER) 

SKIPF(MASTER) 

COPYBF(MASTER,L216) 

RETURN(MASTER) 


' Retrieve the 
■ program from its 
. master tape 


L216. 


EXIT. 

DMP(0, field length) 
— End-of-record 


/ 

17 

Card input 
data 



r Prepare optional 
\ input data files 

{ Execute L216 (DUBFLX) 

r Save optional output 
I data files 


— End-of-file 
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The following list is a typical set of control cards used to execute L216 using the 
relocatable binaries from the program’s master tape. 

Job card 
Account card 


KL,VSN=66XXXX) 

’ r Retrieve the program 

I from its master tape 

’ Prepare optional input files and 
retrieve DYLIB, the DYLOFLEX 
. alternate subroutine library 


LDSET(LIB= DYLIB, 

PRESET=INDEF) 

LOAD(REL216) 

NOGO. 

RETURN(REL216, DYLIB) 

L216. 

* r Save optional output 

• I data files 
EXIT. 

DMP(0, field length) 

— End-of-record 



End-of-file 


fLoad and Execute 

Il216(DUBFLX) 


REQUEST(MASTER, F=I,LB= 
REWIND(MASTER) 
SKIPF(MASTER,2) 
COPYBF(MASTER, REL216) 
RETURN(MASTER) 


6.2 RESOURCE ESTIMATES 

The computer resources utilized (core requirements, tjqies, printed output, time, etc.) 
are a function of problem size. Table 2 shows examples of resources used for some 
problems of varying size. 

6.2.1 FIELD LENGTH 

Program L216 requires a field length of 76 000 octal words to load; execution requires a 
minimum of 132 000 octal words. 

During computation of the inverse of the kernel function matrix, partitions are formed 
whose size depends on the core available in blank common. The number of input/output 
operations varies as the number of partitions. There is a trade-off between central 
processor time and field length for cases involving a large number of aerodynamic 
elements. 
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Table 2.— Examples of Resource Utilization 



Example 

1 

2 

3 

4 

5 

Number of reduced frequencies 

2 

2 

1 

1 

1 

Number of lifting surface boxes 

32 

52 

102 

162 

225 

Number of slender body elements 

6 

10 

0 

0 

0 

Number of modal degrees of freedom 

3 

3 

1 

1 

1 

CPU seconds 

13 

24 

23 

67 

157 

Disk requests 

469 

590 

349 

582 

845 

Disk sectors 

8184 

9150 

7344 

8942 

14836 

Required core 

143000 

149000 

162000 

1 78000 

1 97000 

Printed lines (NPFM = 1) 

1827 

2180 

1697 

2095 

3098 


An equation to determine the recommended field length is: 


Field length = (132000)8+ (120*NTOT)io 


where NTOT is the total number of aerodynamic elements (lifting surface boxes and 
slender body line elements). 

6.2.2 EXECUTION TIME 

The time required to compute unsteady aerodynamic characteristics per specified 
reduced frequency value depends primarily on the number of aerodynamic elements, 
secondarily on the field length. Problem computational time increases roughly as the 
square of the number of elements. Table 2 includes the CPU times for some variously 
sized problems. 

6.2.3 PRINTED OUTPUT 

The number of lines of printed output varies with problem size and specified print 
options. The maximum number of lines printed per reduced frequency is approximately 
10 000 when the checkout print option is specified. Table 2 includes the line count for 
some examples. 
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6.2.4 TAPE DRIVES 


A magnetic tape drive is required if the program is accessed from a master tape. One 
tape drive may be needed for modal input, two for saved aerodynamic and geometric 
data files. 

6.2.5 DISK STORAGE 

L216 uses internal scratch files. Two of these are random access files (TAPES and 
TAPE9); the others are sequential 


6.3 INPUT DATA 

Input data for program L216 chiefly consists of card input data from file INFIL (TAPES) 
optionally supplemented by modal data from file SATAP. 

A summary of the card input data is given at the end of this section. This summary, a 
quick reference for the necessary card input, is included for the convenience of users 
having attained familiarity with the program. 

The alternative forms of unsteady aerodynamic results specified for output have an 
impact on the input card sequence. The order in which cards are input is shown in 
figure 5 . 

6.3.1 FORMAT OF CARD INPUT DATA 

Card data is read in fixed fields and specific columns of the cards as indicated on the 
pages following. The following conventions, with the exceptions noted, are general 
throughout the program. 

1. Floating point variables are read with format ElO.O. 

2. With the exceptions of cards 13.1, 13.2, and 15.1, integer variables are read with 
format 15. 

3. Hollerith variables (keywords, etc.) are read with format AlO. 

Only the first five characters of keywords are checked for recognition. 

6.3.2 CARD INPUT DATA SPECIFICATIONS 

All card sets must be present except as noted. Unspecified columns on any card may be 
used for comments. The underlined five characters in keywords must be present. 


21 
















m 



Figure 5.— (Continued) 
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Figt ire 5. — ( Conchjdec F 
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Card Set 1.0-L216 Data Set Initiation 



KEYWORD/ 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

SDUBLAT 

AlO 

Introduce* L216 card input data aet - must be the first 




card read by L21G 

11-70 



Available for comnents 


The following cards (2.0 through 17.4) are repeated for each data case. 


Card Set 2.0-Title 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

9TITLE 

AlO 

Introduces case title and comments 

11-70 



Available for cocsnents . 




Intended as case title and label for printed output. 
Printed at beginning of output listing for current case 

The number of $TITLE cards 

la unlimited 


Card Set 3.0-Case Identification 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

CASE 

AlO 

Introduces input data for current case. 

11-15 

I CASE 

15 

Case number for current data set 
(1 ^ ICASE < 36) 

16-20 



Not used 

21-30 

CHECKOUT 

AlO 

Requests checkout printing 

This option overides all print flags and a large amount 
of intermediate data is printed. Generally not exercised. 

31-35 

ICOND 

15 

Condition number < = Case number) 

36-70 



Available for comments 




Card Set 4.0-Condition and Configuration Constants 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

FMACH 

ElO.O 

Mach number, tO ^ FMACH < 1.0) 

11-20 

ACAP 

ElO.O 

Reference area (length^) 

Used in calculation of st^^bility derivatives 

21-30 

REFCHD 

ElO.O 

Reference chord (length) 

Used in calculation of section coefficients and st^d^ility 
derivatives 

31-40 

REFSPN 

ElO.O 

Reference semispan (length) 

Used as normalizing factor for generalized force coefficients, 
also appears in expression for rolling moment coefficient. 

If modal input is via SATAP or output aerodynamic data is 

intended for subsequent use in DYLOFLX, value must be unity 
( 8 = 1.0) . 

41-45 

NDELT 

15 

Symmetry flag (Plane Y = 0) 
= 1 symmetric 

“ -1 iuiti-symjnetric 
•= 0 non-symmetric 

46-50 

NP 

15 

Dumber of panels On all lifting (primary and interference) 
surfaces, (0 <_ NP <_ 40) 

51-55 

NB 

15 

Number of slender bodies (0 ^ NB ^ 20) 

If a body has motion in both Y and Z directions (e.g., nacelle). 
It is represented with two geometrically identical slender 
bodies, each to simulate notion in one of the two directions. 

56-60 

NRF 

15 

Number of reduced frequency values for which aerodynamic 
characteristics are required, (0 <_ NRF < 20) 
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Card Set 5.0-Options and Tape Identification 


KEYWORD/ 

VARIABLE 


6-10 

NAIC 

11-15 

NDATA 

16-20 

NSB 


DESCRIPTION 

Flag to aave geometric and aerodynamic data on files 
NGETP and NAETP 

•• 0 Do not aave 

“ 1 Save 


Flag indicating form of modal data 

^ SA array from tape, tabular input, or 
Polynomial (card aet 14) 

^ . Preselected AIC nodes 

(card aet 13) 


Flag to print matrix of dovmwaah factors CD] 
“ 0 Do not print 

“ 1 Print 


Flag to print matrix (cea) of slender body induced downwash 
factors [aPZ], [APY], and body induced downwash. {wg} 

“ 0 Do not print 

“ 1 Print 


Flag to control printout 

• 1 Input geometry, and modal data, and downwash matrix results 

are printed 

• 2 Computed geometry and modal data are printed in 

addition to the data specified above 

• 3 Intermediate results associated with the analysis 

are printed in addition to the data specified 2d30ve 
This option is equivalent to the )ceyword CHECKOUT 
on card set 3 . 


26-30 Not used 

NAIC must be equal to 0 when DUBFLX is used in DYLOFLEX. 









31-40 NAETP AlO File name of the aerodynamic results to be stored. The 

file name must begin in column 31 with an alphabetic char- 
acter and may contain up to seven characters. Each logical 
file contains all the aerodynamic results for one case. The 
first matrix of each logical file contains the control matrix 
for the aerodyneimic data for the current case. 

41-50 NGETP AlO File name of the geometric results to be saved. The file 

name must begin in column 41 with an alphabetic character 
and may contain up to seven characters. Each logical file 
contains all the geometric data for each case. The first 
matrix contains the control matrix for the geometric data for 
the current case. 
































41-45 J5PECS 15 Synnetry flag (riane Z - 0) 

» 0 Mon-aynmetric 

- 1 Symmetric (biplane or jet) 

“ -1 Anti-aymmetric (ground effect) 

46-50 NPC 15 Mode aelector for Alc generation 

■ 0 Alternative 1 - (plunging, pitching, control 

surface and tab rotation) 

• 1 Alternative 2 - (three cambering modes, control 

aurface and tab rotation) 

Omit if NAIC 1 

51-55 NSV 15 Total number of strips on all vertical panels that lie on 

the plane of aymmetry Y •> 0. Data for vertical panels lying 
on the plane Y - 0 must be input before data for other panels. 
56-60 NBV 15 Total number of boxes on vertical panels lying on the plane 

Y = 0. 

61-65 NYAW 15 Flag to control stability derivative calculations 

= 0 Pitch coefficients calculated 

* 1 Yaw coefficients calculated 

If on card set 4 : 


NDELT - 

0: 

NYAW = 

0 or 1 

NDELT - 

1; 

NYAW - 

0 

NDELT = 

-1: 

NYAW - 

1 


NPRl, NPR2 and NPR3 must be 0 if HAIC ^ 1. Strips are numbered consecutively from 
inboard to outboard per panel for all panels in the input order. 






















Card Set 7.0->Gust Data Set 


Omit card set 7.0 if NGUST = 0 (card set 6.0). This option can not be used in DYLOFLEX. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

GUST 

AlO 

Introduces gust specification data 

11-20 

GZRO 

ElO.O 

Gust reference plane dihedral angle (degrees) 

If GZRO » 0.0, gust velocities are in vertical direction 

21-30 

XZRO 

ElO.O 

Gust reference point, location at which sinusoidal gust 
velocity is unity. 

31-40 

VEL 

ElO.O 

Aircraft velocity (length/sec) 

41-50 

WG 

ElO.O 

Vertical gust velocity at XZRO (length/sec) 


Card Set 8.0-Strip-Box Correlation 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

LIM(I,1) 

15 

Number of first box on strip I 

6-10 

LIM(I,2) 

15 

Number of last box on strip I 

11-70 


1215 

Repeat LIM 1 and 2 for all NSTRIP values using as many cards 
as necessary. Each card contains values for seven strips. 


Boxes are numbered consecutively from leading edge to trailing edge per strip< spanwise 
per panel for all panels in the input order. 


Card Set 9.0-Reduced Frequency Data Set 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

REDUCED 

FREQUENCIES 

AlO 

Introduces array of reduced frequencies 

21-70 




Available for cosnnents 
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Card 9.1-Array of Reduced Frequency Values 


COLE. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

RFREQ 

ElO.O 

Reduced frequency k « (radians) 

11-70 


6E10.0 

Repeat RFREQ for all NRF values using as many cards as neces- 
sary. Each card contains up to seven values. 


Geometric arrangements and nomenclature for lifting surfaces and slender bodies are 
illlustrated in figures 6 and 7, respectively. 


The geometry data of cards 10.0 through 12.3 is defined in the reference axis 
system. If used in DYLOFLEX, this reference axis system should be the same as that 
used in INTERP, L215 (ref. 5). Proper aerodynamic modeling may require surfaces to be 
moved up or down, forward or back, or inboard or outboard from their actual location 
(fig. 6). When using the interpolation (SA) arrays from INTERP, control and force point 
locations on a surface, are transformed from the reference axis system to the surface’s 
local axis system. To insure the proper transformation of these points, any shifting of 
the surface must be taken into account using the variables XSHIFT, YSHIFT, ZSHIFT 
on card set 11.0. 

Card Set 10.0-Geometry Data Set 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

GEOMETRY 

AlO 

Introduces geometric input data set. 

11-70 



Available for comments. 


Omit Cards 11.0 through 11.4 if no panels are defined. Repeat all panel data (cards 11.0 
through 11.4) for each panel. All the panels on a lifting surface must be input 
consecutively. The data for the surfaces must be in the following order. 

1. Primary lifting surface panels (e.g., wings, tails, etc.); panels on the Y = plane 
must appear first 

2. Interference lifting surface panels 


32 








Card Set 11.0-Panel Geometric Definition 


COLS. 

KEYWORD/ 

VARIABLE 

1-10 

PANEL 



11-15 

IPN 

16-20 

IDSURF 

21-25 


26-30 

itype 

31-40 

XSHIFT 

41-50 

yshift 

51-60 

ZSHIFT 


DESCRIPTION 


Introduces data for specific panel 


Panel number 

Panels are numbered consecutively as input; Applicable only 
to primary lifting surface panels. 


Logical file number of the SATAP file containing the modal 
data for the surface including this panel. For interference 
surface lifting panels IDSURF “ 0. 


Not used 


Indicates type of lifting surface 
PRIME - primary lifting surface 
INTER - interference lifting surface 


Shift of surface along X axis 


Shift of surface along Y axis 


Shift of surface along Z axis 


Card 11.1 

-Panel Planform Definitioi 

1 



COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 



1-10 

XI 

ElO.O 

X coordinate of panel inboard leading edge corner 



X2 


21-30 I X3 


31-40 


41-50 Y1 
51-60 Y2 


ElO.O 

X coordinate of panel inboard trailing edge corner 

ElO.O 

X coordinate of panel outboard leading edge corner 

ElO.O 

X coordinate of panel outboard trailing edge corner 


Y coordinate of panel inboard edge 

Y coordinate of panel outboard edge 


Card 11.2-Panel Planform Definition (Cont.) 


KEYWORD/ 

COLS . VARIABLE FORMAT 


DESCRIPTION 


21-25 I NC 


26-30 


31-40 


Z coordinate of panel inboard edge 


Z coordinate of panel outboard edge 


Number of chordwise division boundaries 


Number of spanwise division boundaries 


Scale factor for panel deflection modes 
Default; COEFF - 1.0 





















































Card 11.3-Array of Panel Chordwise Division Boundaries 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

TH 

ElO.O 

Chordwise division boundary 

11-70 


6E10.0 

Repeat TH for all NC values using as many cards as 
necessary. Each card contains seven values. 

Fractional values beginning with 0.0 at the panel leading edge, ending with 1.0 at panel 
trailing edge. 


Card 11. 4- Array of Panel Spanwise Division Boundaries 


COLS. 

keyword/ 

VARIABLE 

FORMAT 

DESCRIPTICW 

1-10 

TAU ( I ) 

ElO.O 

Spanwise division boundary 

11-70 


6E10.0 

Repeat TAU for all NS values using as many cards as neces- 
sary. Each card contains seven values. 

Fractional values beginning with 0.0 at the panel inboard edge, ending with 1.0 at 
the panel outboard edge. 


Omit cards 12.0 through 12.3 if no slender bodies are defined. Repeat all slender body 
data (cards 12.0 through 12.3) for each slender body. Slender bodies with Y doublets 
must appear first. 

Card Set 12.0-Slender Body Geometric Definition 


COLS. 

KEYWORD/ 

VARIABLE 

1 

FORMAT 

DESCRIPTION 

1-10 

BODY 

AlO 

Introduces data for specific slender body 

11-15 

IBN 


Body number, (1 < IBN < 20) 

16-20 

IDSURF 

15 

File number of the SATAP file containing the data for this 
slender body. 

21-70 



Available for comments. 
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Card 12.1-Body Axis and Motion Specification 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTIOH 

1-10 

zc 

ElO.O 

Z coordinate of slender body axis 

11-20 

YC 

EIO.O 

Y coordinate of slender body axis 

21-30 

COEFF 

£10.0 

Scale factor for slender body deflection modes 
Default: COEFF -= 1.0 

31-35 

NF 

15 

Number of slender body element endpoints, (2<NF^20) 

36-40 

NZ 

15 

Flag to specify Z doublets (upwash) 
X 0 No Z doublets 
« 1 Z doublets 

41-45 

KY 

15 

Flag to specify Y doublets (sidewash) 
= 0 No Y doublets 

= 1 y doublets 

46-50 

MISBl 

15 

Number of first box on the interference lifting surface 
associated with this slender body 

51-55 

MISB2 

IS 

Number of last box on the interference lifting surface 
associated with this slender body 

If a body has both Z and Y motion (e.g., nacelle), it is represented with two bodies. 

one having Z doublets, one having Y doublets. A single interference lifting surface is 

associated with both bodies 



Card 12.2-Array of Body Element Endpoints 



KEYWORD/ 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

F 

ElO.O 

X coordinates of slender body element endpoints. 

11-70 


6E10.0 

Repeat F for all NF values using as many cards as neoessary. 
Each card may ocntain seven values. 

Ordered from nose 

to tail. 



Card 12.3-Array of Radii at Body Element Endpoints 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 


1-10 

RAD 

ElO.O 

Slender body radii at element endpoints. 
Slender body must be closed, i.e. , RAD(1] ~ 

RAD(NF) = 0.0 

11-70 


EEIO.O 

H^ieat RAD for all NF values using ^ls mar^ cards as 
Each card may ocntain seven values. 

necessary. 




Modal data may be input in any one of three forme: 

1. Coefficients used in polynomial definitions of modal displacements and slopes are 
input on cards. Equations will be written to define the variables. For panels 
(lifting surfaces) the deflections will be normal to the surface, while for bodies the 
deflections will be in the Z- (if any) and^Y- (if any) directions. For the 

panel or body, the deflections in the "NQ” mode are calculated as follows: 

NQ^ “ ^ 2 ARQ(nr.NQ.N.M) (15) 

where t is in a radial direction. The origin of the radius is either at the origin of 
coordinates, when N8 = 0, or at the inboard edge (for the panel) or the axis (for a 
body) if N8 is set to 1 (fig. 8). 

2. Modal deflections, slopes, etc. for each aerodynamic element are input directly on 
cards. 

3. Arrays of modal coefficients are read from input file SATAP and routines within 
L216 implement an appropriate interpolation scheme to determine the required 
modal data at the aerdynamic element control points, defined by geometric data 
previously input. This option requires previous use of the interpolation program 
L215 (INTERP) defined in reference 5, in order to form the modal coefficient arrays 
contained on SATAP. The interpolation routines are defined in the reference. 

Card Set 13.0-Modal Data Set 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

MODES 

AlO 

Introduces modal data input set 

11-70 



Available for comments 
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for panels 



(N8) |^ '+ |z - 


for bodies 



s I 
\ 


\ 

Since each panel is planar, t is a spanwise distance 
in the plane measured from its inboard edge or from 
the origin y = z = 0; t is a radial distance only in 
the sense that each panel may have a different 
dihedral. The use of this radial distance is not 
meant to imply that the panels are curved. 


Figure 8. — Radial Origin Definition 


Omit cards 13.1 through 13.4 if NAIC = 0 (card set 5.0). 
Omit card 13.1 if NPC = 1 (card set 6.0). 

Omit card 13.2 if NPC = 0. 

Card Set 13.1-AIC Panel Mode Selection 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-3 

NOP 

13 

Panel number 

4 

ISl 

11 

= 1 Select plunging mode 

5 

IS2 

11 

= 1 Select pitching mode 

6 

IS3 

11 

= 1 Select control surface rotation 

7 

IS4 

11 

= 1 Select t2d3 rotation 

8 

IS5 

11 

“ 1 Select control surface plunging 

9 

IS6 

11 

1 Select tab plunging 

10 



Not used 

11-70 



Repeat columns 1 to 10 for each panel; i.e. 7 panels 
per card. There are NF (card set 4) panels. 


Card 13.2-AlC Panel Mode Selection (Alternative 2) 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-3 

NOP 

13 

Panel number 

4 

ISl 

11 

Select 1st ceunbering mode («1) 

5 

IS2 

11 

Select 2nd cambering mode (-1) 

6 

IS3 

11 

Select 3rd cambering mode (•!) 

7 

IS4 

11 

Select control surface rotation (>1) 

8 

IS5 

11 

Select tab rotation (“1) 

9 

IS6 

11 

Select control surface plunging (•°1) 

10 

IS7 

11 

Select tab plunging (^1) 

11-70 



Repeat coliimns 1 to 10 for each panel defined, (that is, 
7 panels per card) . There are NP (card set 4 ) panels. 
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Repeat cards 13.3 and 13.4 as pairs for all panels. 

Card 13. 3- Array of Reference Axis Locations at Panel Inboard Edge 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

XHI (1) 

ElO.O 

Elastic axis 

11-20 

XHI (2) 

XIO.O 

Control surface leading edge 

21-30 

XHI (3) 

ElO.O 

Control surface rotation point 

31-40 

XHI (4) 

BIO.O 

Tab leading edge 

41-50 

XHI (5) 

ElO.O 

Tab rotation point 

Fraction values of local chord at the Inboard panel edge 


Card 13. 4- Array of Reference Axis Locations at Panel Outboard Edge 



KEYWORD/ - 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

XHO(l) 

ElO.O 

Elastic axis 

11-20 

XHO{2) 

ElO.O 

Control surface leading edge 

21-30 

XHO(3) 

ElO.O 

Control surface rotation point 

31-40 

XHO(4) 

ElO.O 

Tab leading edge 

41-50 

XHO(5) 

ElO.O 

Tab rotation point 

Fractional values 

of local 

chord at the outboard panel edge 



Omit cards 14.0 through 17.4 if NAIC = 1 (card set 5.0). 


Card Set 14.0-Modal Data Flags and Options 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-5 

NMD 

15 

Total number of modes (1 ^NMD ^ 70) 

6-10 

NTA 

15 

Total number of ARQ values for polynomial mode input 
Omit if NMDIN ^ 0 

11-15 

lAERO 

15 

Flag to save aerodynamic data on file NAETP 
• 0 Do not save 

« 1 Save [d] 'and matrices (xT'is the quasi-inverse) 

“ 2 Save (Q] and fCp] matrices 

- 3 Save [d],' [f] , [q] , [ ^Cp]ma trices 

16-20 

NMDIN 

15 

Type of modal input* 

= 0 Polynomial 

- 1 Integration matrix (B], deflections [h/s], 

first derivative [d (h/s)/d (x/s) ] and second 
derivative [d^ (h/s)/d (x/s) *] from cards. 

« 2 Input SA arrays from file NTPSA 

21-25 

NPFM 

15 

Flag to compute and print pressure coefficients, section 
coefficients, stability derivatives and generalized force 
coefficients 

> 0 Do not compute and print 

- 1 Compute and print 

If lAERO > 1 and NPRM » 0 , the results are computed but not 
printed. 

26-30 

NMTP 

15 

Number of ARQ coefficients for all lifting surface panels 
Omit if NMDIN ^ 0 

31-35 

NMTB 

15 

Number of ARQ coefficients for all slender bodies 
Omit if NMDIN - 0 

36-45 

NTPSA 

AlO 

Name of file containing the SA arrays. File name must begin 
in column 36 with an alphabetic character and may contain 
up to seven characters . Omit if NMDIN / 2 


*h is defined positive in the positive z or y direction (fig. 1) 
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Omit cards 15.0 through 17.4 if NMDIN = 2 (card set 14.0). 

Omit cards 15.0 through 15.2 if NMDIN = 1. 

Omit cards sets 16 through 17.4 if NMDIN — 0. 

Repeat cards 15.0 through 15.2; once for panel data, once for slender body data. 
Card Set 15.0-Polynomial Mode Input Data Set 


COLS. 

KEYWORD/ 

VARIABLE 

FOR.MAT 

DESCRIPTION 

O 

PANEL 

or 

BODY 

AlO 

Indicates polynomial modal data for specific panels, bodies 

11-70 




Available for comments 

Polynomial modal data for panels must precede that for bodies. 


Card 15.1-Polynomial Mode Specifications 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-2 

NR 

12 

Panel/body number 

3-4 

NQ 

12 

Mode number 

5-6 

NA 

12 

Number of ARQ coefficients for mode NQ for panel/body NR. 

7-8 

N8 

12 

Radial origin flag 

>= 0 Use X axis 

• 1 Use panel inboard edge/body axis 

9-10 



Not used 

11-70 


3012 

Repeat the four items in column 1-10 for all ( (NF/NB) *NKD) 
values using as many cards as necessary. Each card may 
contain seven sets of four items. 


Ordering Is required to be ^^Panels (1-NP) Mode(l)| — '^Panels (1-NP) , ModeCNMD) 
or equivalently for slender bodies. The set of numbers must be specified regardless 
of whether or not the panel/body is in motion. 





Card 15.2-Constant Coefficients for Definition of Polynomial Mode Shapes 


Note sign convention shown in figure 1. 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

m 

LARQ(l) 

15 

Exponent of X/S in the expression defining polynomial modes 

6-10 

LARQ(2) 

15 

Exponent of T/S 

11-20 

ARQ 

ElO.O 

Coefficient in the expression defining polynomial modes 

21-60 



Repeat to I NA - three Bets per card - mamber of cards as 
necessary 


If NMDIN = 1, cards 15.0 through 15.2 are replaced by card sets 16.0 through 17.4. 
Repeat cards 16.0 through 16.3 for each mode. 

Card Set 16.0-Tabular Modal Data for Panels 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

i 

1-10 

PANEL MODE 

AlO i 

Introduces tabular modal data set for lifting surface panels 

11-15 

INMD 

15 1 

Mode number, (1 ^ iNMD ^ NMD) 

16-20 

NBOX 

15 

Total number of boxes on all panels 

21-25 

IPLAG 

15 

Flag describing form of input for generalized force integration 
matrix 

* 0 Integration matrix [BQ] input on cards 

- 1 [BQ] calculated internally from normalized box 

c/4 displacements input on cards 


Card 16.1 -Box Integration Matrix 


DESCRIPTION 


ElO.O Generalized force In'tegration matrix data 


6E10.0 Repeat BQ for all NBOX values, using as many cards as necessary. 
Each card contains seven values. 


Form of the data is dependent on IFLAG (card set 16.0) 
IFLAG = 0 Input matrix [BQ] with elements defined by 

IFLAG - 1 Input matrix [ ] , c/4 displacements 



















































Card 16.2-Box Displacements 


Note sign convention shown in figure 1. 



KEYWORD/ 



COLS. 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

H 

ElU.O 

Normalized box 3c/4 displacements J 

11-70 


6E10.0 

Repeat H for all NBOX values using as many cards as necessary. 



■ 

Each card contains seven values. 

Values for boxes on interference lifting surfaces may be input as zero, since interference 
surfaces carry no lift due to their own motion. 


Card 16.3-Derivatives of Box Displacements 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

DH 

ElO.O 

Derivative of normalized box 3c/4 displacement 




d (h3c/4/s) /d (x/s) 

11-70 


6E10.0 

Repeat OH for all NBOX values using as many cards as necessary. 




Each card contains seven values. 


Values for boxes on interference surfaces may be input as zero. 


Card Set 17.0-Tabular Modal Data for Bodies 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

BODY MODE 

AlO 

Introduces tabular modal data set for slender bodies 

11-15 

INMD 

IS 

Mode number, (1 < INMD < HMD) 

16-20 

NBE 

-J 

Total number of line elements on all slender bodies 


Card 17.1 -Body Element Integration Matrix 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

BQ 

ElO.O 

Generalized force Integration matrix data 

11-70 


6E10.0 

Repeat BQ for all NBE values using as many cards as necessary. 
Eacn card contains seven values. 


Form of data is dependent on IFLAG (card set 16.0) 

IFLAG = 0 Input matrix [BQ] with elements defined by 


BQ^i - 1, [r.^ 

IFLAG = 1 Input matrix J 
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Card 17.2-Line Element Displacement 


Note sign convention shown in figure 1. 


COLS. 

KEYWORD/ 

VARIABLE 

1 

FORMAT 

DESCRIPTION 

B 

H 

ElO. 0 

Normalized line element midchord displacement j 

11-70 


6E10. 0 

Repeat H for all NBE values using as many cards as necessary. 
Each contains seven values. 


Card 17. 3-First Derivative of Line Element Displacement 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-10 

DUl 

ElO.O 

Derivative of normalized line element midchord displacement 

11-70 


6E10.0 

Repeat DHl for all NBE values using as many cards as necessary. 
Each contains seven values- 


Card 17.4-Second Derivative of Line Element Displacement 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

■ 

DH2 

ElO.O 

Second derivative of normalized line element midchord dis- 
placement d* 

■ 


6E10. 0 

Repeat DH2 for all NBE values using as luany cards as necessary. 
Each contains seven values. 


Card Set 18.0-L216 Data Set Termination 


COLS. 

KEYWORD/ 

VARIABLE 

FORMAT 

DESCRIPTION 

1-20 

$QUIT 

DUBLAT 

DATASET 

AlO, lOX 

Terminates L216 card input data. 
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Requirements 
or Function 


Key Words and/or Variables 


Card Format 


Reference 
Card Set (CS) 



$DUBL 

A10 

1.0 


$TITL 

AlO 

2.0 


CASE CHECK 

ICASE ICOND 

A10.I5, 
AlO, 15 

3.0 

Condition, 

Configuration 

Specs. 

FHACH ACAP REFCND REFSPN NDELT NB 

HP NRF 

4E10. 0,415 

4.0 


Options 


NDSV NDATA NDBG 
NAIC NSB 


NAETP NGETP 


5I5,5X,2A10 


5.0 


AIC Specs., 


NSTRIP NBL 


KPR2 


JSPECS NSV 


NYAW 


6IS,5X,6I5 


Options NBF NPRl NPR3 NGUST NPC NBV 

Gust; NGUST-1 GUST GZRD XZRO VEL WG j A10.4E10.0 7.0 

Strip-Box LIH(I.I) , Repeat to I = NSTRIP \ 1415 8.0 


Correlation 

LIM(I,2) 




REDUC 

AlO 

9.0 

k values 

RFREQ , Repeat to NRF 

7E10.0 

9.1 


GEOHE 

AlO 

10.0 

Panel Specs. 
HP )« 0 

PANEL IPN XSHIFT YSHIFT ZSHIFT 

IDSURF ITYPE 

Repeat 

to 

HP 

AlO, 215, 5X, 
A5.3E10.0 

11.0 

Panel Planform 
Definition 

XI X2 X3 X4 Y1 Y2 

6E10.0 

11.1 

Z1 11 NC COEFF 

NS 

2E10.0.215 

ElO.O 

11.2 

Chordwlse Divisions 

TH . Repeat to NC 

7E10.0 

11.3 

Spanwise Divisions 

TAU , Repeat to NS 

7E10.0 

11.4 




















































NWDIII - 1 NWIN - 0 » X l«rC 


Requirements 
or Function 

Key Words and/or Variables 

Slender Body Specs. 
NB 0 

BODY I BN 

IDSURF 

Body axis, motion 
Specs . 

ZC YC COEFF NF NY M1SB2 

NZ MISBl 

Body Element 
Endpoints 

F , Repeat to NF 

Body Radii at 
Endpoints 

RAD , Repeat to NF 

j 


AIC Specs. I 
HPC ^ 0 

f AIC Specs. II 
NPC = 1 

Ref. Axis.inbd. 

Ref. Axis.Outbd. 

odal Data Specs, 
nd Options 


Poly. Mode Specs. 
Coefficients 

Box Tabular Modal 
Specs. & Options 

I 

Box Integ. Matrix 

Box H3c/4/s 

Box dH3c/4/s/dx/s 

Body Tabular 
Modal Specs. 

Line Integ. Matrix 

Line Hc/2/s 

Line dHc/2/s/dx/s 

Line d*Hc/2/s/dx/s' 


[N0P.IS1.IS2,IS3,IS4,IS5,IS6], Repeat to NP 


[N0P,IS1,IS2,1S3,IS4.IS5,IS6,IS7], Repeat to NP 



AlO 


7(13,611) 


7(13,711) 



NTA NMD IN NHTP NTPSA 


PANEL/BODY 


[NR,NQ,NA,N8,-], Repeat to NP/NB 


LARQl ARQ , Repeat to NA 
LARQ2 


PANEL I NMD NFLG 
NBOX 


BQ , Repeat to NBOX 


H , Repeat to NBOX 


DH , Repeat to NBOX 


BODY I NMD 


AlO 

412, 2X, 3012 
2I5,E10.0 


























6.3.4 MAGNETIC FILES INPUT DATA 


Program L216 utilizes one magnetic file (disk or tape) for modal data input. 
Interpolation arrays composed of modal deflections and slopes are input on file SATAP. 
This file must have been created after previous execution of the program described in 
reference 5. Order and format of the matrices on file SATAP are shown in figure 9. 


MATRIX SIZE 
(NNODES X NTMODE) 
(NNODES X NTMODE) 
(NNODES X NTMODE) 
(NNODES X NTMODE) 
(NNODES X NTMODE) 
(NNODES X NTMODE) 
(NNODES X 6) 

See formula in 
reference 5. 



EOP 

EOF 


i^*^ surface in i^^ file 


Optional SA 
EOF 


Note: 

6 , d , 6 are in degrees 
A y ^ 

NNODES = number of nodes 
NTMODE = number of modes 


Figure 9. —Map of Magnetic File 54 TAP 


6.4 OUTPUT DATA 

Output data from program L216 consists of printed results, and geometric and 
aerodynamic data saved on magnetic files NGETP and NAETP at the user’s request. 
Saved files are used subsequently in DYLOFLEX. 
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6.4.1 PRINTED OUTPUT DATA 


The specific data and the extent of printed output are dependent on user specified 
options. The major blocks of printed output are described below. 

1. Echo print of input options, flags, title, file names, etc 

2. Geometric data input for lifting surfaces 

3. Geometric data input for slender bodies 

4. Modal data 

• Normalwash due to primary lifting surface motion 

• Normalwash induced by slender bodies 

• Augmented normalwash 

5. Aerodynamic characteristics/per mode 

• Pressure difference coefficients 

• Section normal force and moment coefficients 

• Slender body forces 

• Stability derivatives 

• Generalized airforce coefficients 

6. Intermediate geometric and aerodynamic data (optional) 

6.4.2 MAGNETIC FILES OUTPUT DATA 

Selected geometric and aerodynamic data can be saved at the user’s request on files 
"NGETP” and "NAETP”, respectively, for subsequent use in DYLOFLEX. Except for 
the quasi-inverse matrix (D ^), the files are writen in READTPAVRTETP format^; an 
end-of-file terminates the matrices stored for a case. The quasi-inverse matrix is written 
with binary READ/WRITE format in multiple records terminated by a one word record 
containing (lOHQUASI-END.). 


^Clemmons, R.E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface with FLEXSTAB. NASA contract NAS 1-13918, BCS-G0701, September 1975. 
(internal document) 
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Saved Geometry Tape NGETP 


The format of all geometric data matrices on file NGETP is described in figure 10. 
Contents of the first matrices on file are shown in figures 11 through 15. The 
interference surface-slender body correlation matrix is described in figure 16. 


SIZE 

15 

8*NTB 

Nl 

N2 

N3 

NBOX 

each 


NBOX 

each 

NTOT 

NTOT 

NTOT 

NBE 

NBE 

NBE 

NBE 


MATRIX 

NO. DESCRIPTION 


1 

MPC 

2 

BODTAB 

3 

MPSG 

4 

MISG 

5 

MSBG 

6,7,8 

EV 


PV 


ZV 

9,10,11 

X 


Y 


ZZ 

12 

MBA 

13 

MDA 

14 

MSBC 

15 

SDELX 

16 

X 

17 

RO 

18 

ROP 


EOF 


Geometry control matrix (See Figure 11) 

Body table matrix (See Figure 12) 

Primary lifting surface geometry matrix 
(See Figure 1 3 ) 

Interference lifting surface geometry 
matrix (See Figure 14) 

Slender body geometry matrix 
(See Figure 15) 

Arrays of X,Y, and Z coordinates of 
the 1/4 chord points of all boxes 

Arrays of X,Y, and Z coordinates of 
the 3/4 chord points of all boxes 

Aerodynamic element area 

Dihedral angles of aerodynamic elements 

Interference surface box-slender body 
line element correlation (See Figure 16) 

Body element lengths 

Body element midpoint X location 

Body element midpoint radius 

Slope of midpoint radius 


Figure 10.— Geometry Definition Tape "NGETP" 
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SIZE 


15 


NO. 

1 

2 

3 

a 

5 

6 

7 

8 

9 

10 
1 1 

12 
1 3 

1*4 

15 


ITEM 

ICASE 

ICOND 

IGEOM 

REFCHD 

REFSPN 

ACAP 

FMACH 

NTBB 

NIBB 

NBE 

NTB 

NB 

NPLS 

NPIS 


DESCR IPTION 

Case Number 
Condition Number 

Geometric data successfully saved (=1) 

Reference chord 
Reference semispan 
Reference Area 
Mach Number 

Number of boxes on thin bodies (primary lifting 
surfaces) 

Number of boxes on interference surfaces 
Number of slender body elements 
Number of panels on thin bodies (primary 
lifting surfaces) 

Number of slender bodies 

Number of thin bodies (primary lifting 

surfaces) 

Number of interference lifting surfaces 
Not Used 


Figure 1 1. —Geometry Control Matrix [MFC] 
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SIZE = 8*NTB 


The matrix contains one column for each body I. Elements in a column 
are specified below. 

element^ fiESCRIPTIQIJ 


(1,1) 

ITYPE 

= 1, Slender Body 
= 3, Thin Body 

(2,1) 

IPOS 

= 0, Off plane of symmetry 
= 1, on plane of symmetry 

(3,1) 

Xo 

Origin of local axis system 

(4,1) 

Yo 

for this body (Reference 

(5,1) 

Zo 

axis coordinate) 

(6,1) 


Dihedral angle for this body 

(7,1) 

IFIRST 

Pointers to this body's first and 

(8,1) 

LAST 

last boxes/segments in the total array 
of structural boxes/segments 


NOTE; A thin body corresponds to a lifting surface panel. 


Figure 12. —Body Table Matrix [BODTABJ 
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NLIFS 

Size= [4 + 5*NLiFS+ (NSUBL (2,l)*9) ] 

I = 1 


No. 


Item 


Description 


1 

ICASE 

2 

ICOND 

3 

NLIFS 

4 

NBXL 

5 

NSUBLdJ) 

6 

NSUBL(2,1) 

7 

NSUBL(3,1) 

8 

NSUBL(4,1) 

9 

NSUBLf 5,1) 

10 

IPN 

11 

XI 

12 

X2 

13 

X3 

14 

X4 

15 

Y1 

16 

Y2 

17 

Z1 

18 

Z2 


; i 

' 



■■ 


Case number 

Condition number 

Number of primary lifting surfaces 

Number of boxes on all primary lifting surfaces 

File number of SA array for first surface 

Number of panels on first surface 

First box number on first surface 

Last box number on first surface 

Flat = 1 if surface on Y = 0 plane 

Panel number 

^ Panel coordinates (see fig. 5) 


Elements from 10 through 18 are repeated for 
each of the NSUBL(2,1) panels 

Elements 5 through (9*NSUBL(2,1 ) + 9) are 
repeated for all NLIFS surfaces 


Figure 13. —Primary Lifting Surface Geometry Matrix [MPSG] 
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MINTS 

Size= [4 + 5*NINTS+ ^ (NSUBL (2,l)*9)] 
1=1 


No. Item 


1 

ICASE 

2 

ICOND 

3 

MINTS 

4 

NBXL 

5 

NSUBL (1,1) 

6 

NSUBL (2,1) 

7 

NSUBL (3,1) 

8 

NSUBL (4,1) 

9 

NSUBL (5,1) 

10 

IPN 

11 

XI 

12 

X2 

13 

X3 

14 

X4 

15 

Y1 

16 

Y2 

17 

Z1 

18 

Z2 

/ 

4 

^ / 


✓ 

/ 

/ 


Description 

Case number 
Condition number 

Number of interference lifting surfaces 
Number of boxes on all interference lifting surfaces 
File number of SA array for first surface 
Number of panels on first surface 
First box number on first surface 
Last box number on first surface 
Flag = 1 if surface is on Y = 0 plane 
Panel number 
\ 


y Panel coordinates (see fig. 5) 


Elements 10 through 18 are repeated for each of 
the NSUBL (2,1) panels 

Elements 5 through (9 + 9*NSUBL (2,1)) are repeated 
for all MINTS surfaces 


Figure 14. —Interference Lifting Surface Geometry Matrix [MiSG] 



Saved Aerodynamics Tape NAETP 

Contents of the aerodynamic data file NAETP are dependent on the user specification of 
lAERO (card set 14.0). The number of matrices written on NAETP in each instance is 
shown in table 3. The formats of these matrices for the cases when lAERO = 1 or 
lAERO = 2 are the shown in figures 17 through 19. If lAERO == 3, the result is 
equivalent to requesting both lAERO = 1 and lAERO = 2. 

Table 3.~Number of Matrices Stored on NAETP 


lAERO 

Number of matrices 

0 

None 

1 

(2*NBOX + N + 4} NRF 

2 

(NMD + 2) NRF 

3 

(2NBOX + N + NMD + 6) NRF 


Note: For lAERO = 3, the matrices for lAERO = 1 are followed 

by the matrices for lAERO = 2 for each reduced frequency. 
The matrix sizes and contents are described in figures 16 
through 18. 


6.5 RESTRICTIONS 

Restrictions on the theoretical development and on idealization of configurations have 
been previously discussed in section 4.2. Restrictions on program variable sizes have 
been defined in section 6.3. The aerodynamic influence coefficient formulation (AIC) is 
obtained for strips on lifting surfaces only-slender body elements have been ignored. 
The AIC option in this program does not produce results that are suitable for 
subsequent use in DYLOFLEX. 


6.6 DIAGNOSTICS 

Errors detected by L216 generate diagnositc messages. Fatal errors cause the current 
case to be terminated and a new case to be initiated. Thirty-three fatal error messages 
exist in the program coding. The first line of a fatal error reads: 

***** FATAL ERROR** 

and the last line is: 

CURRENT CASE WILL BE TERMINATED 
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lAERO = 1 


NO. 

1 

2 

(NBOX+1) 

(NBOX+2) 

CNBOX+2+N) 

(NBOX+3+N) 

(2*NB0X+3+N) 

(2*NBOX+4+N) 


SIZE 


DESCRIPTION 


(14+NKVAL) 

(2*NBE) 


(2*NBE) 

VARIABLE 


VARIABLE 


( 2*NBOX+l) 


3 


1 


MAC 

[F], 

• 

• 

• 

[f] 

'• ■'nbox 


• 

• 



• 

• 

• 

^^^tInbox 

DEND 




Aerodynamic Control Matrix 
See Figure 19 

1st row of the induced 
downwash factors matrix 


Last row of the induced 
downwash factors matrix 

1st record of the L matrix 
portion of the quasi-inverse 
matrix 


Last record of the L matrix 
portion of the quasi-inverse 
matrix 

1st record of the trapezoidal 
section of the quasi-inverse 
matrix 


Last record of the trapezoidal 
section of the quasi-inverse 
matrix 

Record to terminate the 
quasi-inverse matrix. 

DEND = (lOHQUASI-END. ) 

Matrices 2 through (2*NBOX+N+4) 
are repeated for all NKVAL 
reduced frequencies 

EOF 


Figure 17.— Aerodynamics Results Matrix (IAER0=1 ) 
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I AERO 


2 


NO. SIZE 

1 (14+NKVAL) 

2 2*(NBOX+NBE) 

NMD+ 1 2 * ( NBOX+NBE ) 

NMD+2 2* NMD* NMD 


MAC 


{ACp}i 


{ACp} 


NMD 


[Q] 


DESCRIPTION 


Aerodynamic data control 
matrix (see Fig. 19) 


Column of pressure 
coefficients for the 
first mode 

Column of pressure 
coefficients for the 
NMD mode 

Generalized airforces [q] 

Matrices 1 through (NMD+2) 
are repeated for all NKVAL 
reduced frequencies 


EOF 


Figure 18.— Aerodynamic Results Matrix (IAER0=2) 
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SIZE = 14 ♦ Number of reduced frequencies 

NO. ITEM DESCRIPTION 

Case Number 
Condition number 
Flag for saving aerodynamic data 
Mach number 

Number of reduced frequencies 
Number of lifting surface boxes 
Number of slender body elements 
Number of aerodynamic modes 
Reference chord 
Reference semi-span 
Reference area 
Reference length 
Not used 
Not used 

1st reduced frequency 

• 

Last reduced Frequency 

Figure 1 9.— Aerodynamic Results Control Matrix 
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6.6.1 L216 ERROR CODES 


Code Example message 

1 PREMATURE END-OF-FILE FOUND ON FILE nnn 
EXECUTION TERMINATED 

2 UNRECOGNIZED KEYWORD aaaaa IMPROPER INITIALIZATION OR 
TERMINATION 

Valid 5-character keywords are $DUBL, CASE and $QUIT. See card 
input description (sec. 6.3) for proper intialization or termination of 
data sets. 

3 TOO MANY CASES DEFINED MAX ALLOWED IS 36 
EXECUTION TERMINATED 

4 FATAL ERROR FOUND IN OVERLAY nnn 
EXECUTION TERMINATED FOR CASE nn 
FIND NEXT CASE 

5 INVALID NUMBER OF REDUCED FREQUENCIES nnnnnn 

6 INVALID REFERENCE CHORD nnnnnnnnnnnnnnn 

7 INVALID REFERENCE SEMISPAN nnnnnnnnnnnnn 

8 INVALID REFERENCE AREA nnnnnnnnnnnnn 

9 NO. PANEL OR BODY DEFINED - CHECK DATA 

10 NO. OF PANELS nnnn EXCEEDS MAXIMUM (40) ALLOWED 

11 NO. OF BODIES nnnn EXCEEDS MAXIMUM (20) ALLOWED 

12 NO. OF STRIPS nnnn EXCEEDS MAXIMUM (70) ALLOWED 

13 INVALID SYMMETRY OPTION NDELT = nnnnnn 

14 INVALID MACH NUMBER (MIN = 0.0, MAX = 0.9999) 

15 INVALID YAW MOTION FLAG NYAW = nnn 

16 KEYWORD MISSING FOR MODAL INPUT DATA 

17 INVALID X-COORDINATES FOR PANEL nnnnn 

18 SPANWIDE OR CHORDWISE BOUNDARIES EXCEED MAXIMUM (30) 
ALLOWED 
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19 NUMBER OF BOXES DEFINED EXCEEDS 400 

20 NUMBER OF BODY DIVISIONS EXCEEDS MAX (20) ALLOWED 

21 INVALID BODY ENDPOINTS FOR BODY nnnnn 

22 INVALID CAMBER MODE AXIS DEFINED FOR AIC MODE nnnnn FOR 
STRIP nnnnn 

23 WRONG NUMBER OF BOXES SPECIFIED nnnnn 
PROGRAM FOUND nnnnn 

24 WRONG NUMBER OF BODY ELEMENTS nnnnn 
PROGRAM FOUND nnnnn 

25 FETAD ERROR CODE = nnnnn FOUND WHILE DEFINING SA ARRAY 
TAPE nnnnnnn ERROR CODE INTERPRETATION 

= 1 ILLEGAL TAPE NAME OR NUMBER 

= 2 BUFFER TOO SMALL 

= 3 TOO MANY FILES (MAX 49) 

26 FETAD ERROF CODE = nnnnn FOUND WHILE DEFINING NTPQ 
TAPE aaaaaa 

27 READTP ERROR IRR = nnnnn WHILE READING 1ST MATRIX OF SA 
ARRAY 

28 NUMBER OF MODES REQUESTED nnnnn 
NUMBER OF MODES IN SA ARRAY nnnnn 

29 ERROR FOUND IN COMPUTING MODAL DATA FROM SA ARRAY 
AINTG ERROR CODE WORD IS * 

♦ABCDEFGXX MEANS ERROR NO. XX IN 

ROUTINE ABCDEFG 

♦AINTG 1 MEANS INTERPOLATION COEFFICIENT 

ARRAY TYPE UNRECOGNIZABLE 

30 WRTETP ERROR ENCOUNTERED IN WRITING THE nnnTH MATRIX 
ON FILE aaaaaaa. 

ERROR CODE IRR = nnnnn 

31 PANELS ON SURFACE nnn WERE NOT INPUT CONSECUTIVELY 
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32 


FLAG WAS SET TO SAVE AERODYNAMIC RESULTS INVALID TAPE 
NUMBER (NAETP) = aaaaaa WAS GIVEN 


33 

6.6.2 

Code 

0 

1000 

2 

3 

3000 

4 

5 

6 

7 

6.6.3 

Code 

0 

1000 

2 


INVALID FILE NAME 
READTP ERROR CODES 

Description 


No errors are detected during reading. 

+ I An FSF error occurred, where I is the number of file marks remaining l(j 
be skipped when an end-of-information was encountered. 

The number of matrices or files to be skipped before reading starts is less 
than zero. 

The dimensioned number of rows in the matrix is less than zero. 

+ 1 An FSR error occurred, where I is the number of records remaining to be 
skipped when either an end-of-file or end-of-information was 
encountered. 

Number of rows in the matrix is greater than the dimensioned row size 
in the program. 

The name check failed. 

The number of rows in the matrix (M) times the number of columns (N) 
is greater than the buffer size, or M * N^ 0. 

An end-of-file was read. If it occurs while reading the matrix ID. no 
information is stored in the user’s area. If it occurs while reading the 
matrix, the ID information will be stored. Note that the records will 
always be in pairs, and an end-of-file should always be encountered with 
the ID record. 

WRTETP ERROR CODES 

Description 

No errors are detected during writing 


I An FSF error occured, where I is the number of file marks remaining to 
be skipped when an end-or-information was encountered. 

The number of matrices or files to be skipped, before writing starts, is 
less than zero. 
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3 

3000 + I 
6 


The dimensioned number of rows in the matrix is less than or equal to 
zero. 

The actual number of rows is greater than the dimensioned number of 
rows in the matrix. 

The number of rows in the matrix (M) times the number of columns (N) 
is greater than the buffer size. 
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7.0 SAMPLE PROBLEM 


The configuration chosen for the example problem is illustrated in figure 20; included 
are a wing, strut, nacelle, and fuselage. There are three modes of motion: plunge, pitch, 
and roll. The wing is represented by two panels, one having nonzero dihedral. The 
nacelle is idealized as two slender bodies to accomodate both Z and Y motions. Each half 
of the interference shells is represented with three panels. The polynomial modal input 
option is utilized. The idealization is entirely demonstrational and is not to be regarded 
as an optimum one. 


7.1 INPUT DATA 

7.1.1 PARAMETERS AND GEOMETRY 

k = 0, 0.5 A = 6.4 

M = 0.85 Cr - 1.5 

8 = 1.0 s = 1.0 

e = 0 

Nacelle, Slender Body 1 


Yc = 2.0, Zc = -0.5, 

Y doublets 

Element endpoints 

2.0, 2.5, 3.25, 4.5 

Radii at endpoints 

0.0, 0.5, 0.5, 0.0 

Nacelle, 

Slender Body 2 

Yc = 2.0, Zc = 0.5, 

Z doublets 

Element endpoints 

2.0, 2.5, 3.25, 4.5 

Radii at endpoints 

0.0, 0.5, 0.5, 0.0 

Fuselage, Slender Body 3 

Yc = 0.0, Zc = 0.0, 

Z doublets 

Element endpoints 

0.0, 2.0, 4.0, 6.0 

Radii at endpoints 

0.0, 1.0, 1.0, 0.0 


The panel geometry is summarized in table 4. 





Table 4. -Sample Problem Panel Geometry 


Panel 


><1 


^3 

— 

V 


^2 

Zl 

^2 

Chord wise 

Spanwrse 

Wing 

I 

2.0 

4.0 

2.5 

4.0 

0.866 

2.0 

0.5 

0.5 

0.0.0.5.1.0 

0.0,0.5. 1.0 


2 

2.5 

4.0 

3.0 

4.0 

2.0 

3.0 

0.5 

1.0 

0.0, 0.5, 1,0 

0,0. 0.5, 1.0 

Strut 

3 

2.5 

4.0 

2.5 

4.0 

2.0 

2.0 

0.5 

0.0 

O.D,0.5J,0 

0.0, 1.0 


4 

2.5 

4.0 

2.5 

4.0 

2.0 

2.443 

0.0 

-0.25 

0.0,0.5,1.0 

0.0. 1.0 

Nacelle 

5 

2.5 

4.0 

2.5 

4.0 

2.433 

2.433 

-0.25 

-0.75 

0.0, 0.5, 1.0 

0.0, 1.0 

interference 

6 

2.5 

4.0 

2.5 

4.0 

2.433 

2.0 

-0.75 

-I.O 

0.0.0.5,1.0 

0.0, 1.0 

surface 

7 

2.5 

4.0 

2.5 

4.0 

2.0 

1.567 

-1.0 

-0.75 

0.0, 0.5, 1.0 

0.0, 1.0 


8 

2.5 

4.0 

2.5 

4.0 

1.567 

1.567 

-0.75 

-0.25 

0.0, 0.5, 1.0 

0.0, 1.0 


9 

2.5 

4.0 

2.5 

4.0 

1.567 

2.0 

-0.25 

0.0 

0.0, 0.5, 1.0 

0.0, 1.0 

Fuselage 

. 

10 

2.0 

4.0 

2.0 

4.0 

0.0 

0.866 

1.0 

0.5 

0.0, 0.5, 1.0 

0.0, 0.5, 1.0 

interference 

n 

2.0 

4.0 

2.0 

4.0 

0.866 

0.866 

0.5 

-0.5 

0.0, 0.5, 1.0 

0.0, 0.5, 1.0 

surface 

12 

2.0 

4.0 

2.0 

4.0 

0.866 

0.0 

-0.5 

-1.0 

0.0, 0.5. 1.0 

0.0,0.5,1.0 
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7.1.2 MODAL DATA 


Plunge, Mode 1 


Panels^ 

ARQnR, 1,0,0 


Bodies 

Ai^Ql,l,N,M 

= 0.0 


ARQ2, 1,0,0 

= ARQ3^ 1,0.0 =10 



Pitch, Mode 2 (About x = 0) 

Panels^ 

ARQnr,2,i,o 


Bodies 

ARQi,2,N,M 

= 0.0 


AHQ2,2,1,0 

= ARQ3 ,2,1,0, - 1-0 



Boll, Mode 3 

Panels^ 

ARQnr,3,o,o 

= ^>nr'“'^nr +2Jnr"‘'*-^nk 


ARQnR,3,0.1 

= 1.0, N8 = 1 

Bodies 

ARQi^,o,o 

= 0.5 


ARQ2,3,0,0 

= 2.0 


ABQ3,3,N,M 

= 0.0 


^Except NR = 3, 5, 8, 11 


2^ - r , -Zl 

y - tan _ Yj 


‘NR = 1 to 9 


Boeing Commercial Airplane Company 
P.O, Box 3707 

Seattle, Washington 98124 
May 1977 
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duQl ftr 
I l^Lt 
nut 
nut 

I iriE 

n rit 
nut 
nut 
nut 
M u E 
nut 
nut 
nut 

mi f 
nut 
nut 

I I I Lt 
f UvE 
HI UE 
inuE 
H i TlE 
HITlE 
CAiE 


PPOCBftH L?16 S*«PlE PP061EM WING, ilSvJt, NACELLE. FuSELA&E 

Twelve panels I 2-W|MC ,1-StPUT ,6-NAC£U£ INTEPf EPENCE, 

J-E LISE lace in IEKK PENCE! 

ruPEE SLEnOEP auoits Ev-NACELLE , /-nacelle, fuselage l 

seventeen STAIPS ( S-U (NC , l-JTfluT .C-NACELLE [NTfflf EflgNCE , 
i- F use LACE INTEAF f PENCE I 

THjPTr fnufi BniEi le-NiNc, e-stput , iz-nacelle l Nr£«f erenCe, 
J^-FUSElAt£ jNfERf EPENCEJ 

NINE 1 INC ElENCnTS O-EICN StCNOEP flOorl 

THREE OOF I PLUNCf , PI TCH.POLL/ 

PITCH REF X-0, ROLL PEF Y ■ 0, / • 0 

POLYNOMIAL MODE INPUT 


I 


o.aso 


6. A 00 

1.50 

1 

0 

0 

0 2 

l 7 



0 0 

1 

2 

3 

A 5 

1 s 

16 

17 la 19 

29 

30 

31 32 33 

REOOCEO 

PEQUEfiCIES 


0.00 


0,50 


CtO«ElRV 




panel 


1 

( 

2.00 


A .00 

2.50 

0.50 


0.50 

3 

0.00 


0.50 

1.00 

0.00 


0.50 

1 .00 

panel 


2 

f 

2. SO 


A .00 

3.00 

0. SO 


1.00 

3 

0.00 


0.53 

1.00 

O.QO 


0,50 

1 .00 

PANEL 


1 

( 

2.50 


A. 00 

2.50 

0. SO 


0.00 

3 

0.00 


0.50 

1.00 

O.OO 


1.00 


panel 




2. SO 


A. 00 

2.50 

0.00 


-0.250 

3 

0.00 


0.50 

1.00 

0.00 


1 .00 


PANEL 




2.50 


A. 00 

2.50 

-0.2S 


-0,750 

3 

0.00 


0,50 

1.00 


I 

1.00 

SAVAEP 


iO 

JA 




\ 

SAVGEO 
0 0 

a q 

22 23 


12 


0 

11 

2S 


2.00 


0 

IS 

27 


prime 


1.00 


V.OO 2. 

l.OO 


INTER 


A. 00 2.00 

1.00 


INTER 


2.00 


2.A33 


2.A JJ 


A .0 
S .0 
6.0 
a.o 

B .0 
0.0 
9.0 
9-1 
10 .0 

I ( .0 
1 1 . 1 
11.2 
M .S 

I I .A 

I I .0 
11.1 
11.2 
U .S 
n . A 

II .0 
1 1 . 1 
U .2 
n . 3 
n .A 

I L .0 

n. 1 

I I .2 

11.3 

I i .A 
U-0 

II . 1 
11.2 
11.3 


1.00 


7.2 SAMPLE PROBLEM INPUT AND OUTPUT 



0.00 

PANEL 

1.50 
-0.75 

□ .00 
0.00 
PArjEL 
2. 50 
-l.OO 
0.00 
0.00 
PANEL 

2.50 
-0.7 5 

0.00 
0 .00 
PANEL 
2.50 
-0.25 
0.00 
0.00 
PANEL 
2.00 
1.00 
0.00 
0.00 
PANEL 
2.00 
0.50 
0.00 
0.00 
panel 
2.00 
-0.50 
0.00 
0.00 
eoor 


,00 

,00 


0 . 

0. 

0 . 

5. 

- 0 . 

0 . 

0. 

5. 

-1. 

0. 

0. 


.00 

.50 

.00 

.00 

.50 

.50 

.50 

00 

50 

50 

50 

00 

UO 

50 

50 


INIER 

2.50 <>.00 

3 2 

1.00 

INTER 

2.50 5.00 

3 2 

t.OO 

INTER 

2.50 5.00 

3 2 

1.00 

INTER 

2.50 5.00 

3 2 

1.00 


2.533 


1.00 


INTER 


2.00 

3 

1.00 

1.00 

1 

2.00 

3 

1.00 

1.00 

1 

2.00 

3 

1.00 

1.00 


1.00 


5.00 0.866 

1.00 


0.000 


y NOTION nacelle 


11.5 
11 .0 
11.1 
11.2 
11.3 
11.5 
11.0 
11.1 
11.2 
11.3 
11.5 
11.0 
11.1 
1 1.2 
11.3 
11.5 
11.0 
11.1 
11.2 
11.3 
11.5 
11 .0 
11. 1 
11.2 
11.3 
11.5 
1 1.0 
11.1 
11.2 
11.3 
11.5 
11.0 
11.1 
11.2 
11.3 
11.5 
12.0 


-0.500 

2.00 


1.00 5 

0 

1 

11 


22 




12.1 

2.0 

2.50 

3.25 

5.50 









12.2 

0.00 

0.50 

0.50 

0.00 









12.3 

eoov 

2 

1 

NUriUN NACELLE 








12.0 

-0.500 

2.00 


1.00 5 

1 

0 

11 


22 




12.1 

2.0 

2.50 

3.25 

5.50 









12. 2 

0.00 

0,50 

0.50 

0.00 









12.3 

BOOT 

3 

FUSELAGE 









12.0 

0.00 

0.00 


1.00 5 

1 

0 

23 


35 




12.1 

0.00 

2.00 

5.00 

L.OO 









12.2 

0.00 

1.00 

1.00 

0.00 









12.3 

NODES 

PULYNONIAL 

COEFFICIENTS 









13.0 

3 

50 3 

0 1 

35 6 









15.0 

PANEL 












15.0 

1 1 1 

2 1 1 

3 1 0 

5 1 1 

5 

1 0 


6 

1 1 

7 

1 

1 

15.1 

8 1 0 

9 1 1 

10 1 1 

11 1 0 

12 

1 1 


1 

2 1 

2 

2 

1 

15.1 

3 2 0 

5 2 1 

5 2 0 

6 2 1 

7 

2 1 


B 

2 0 

9 

2 

1 

15.1 

10 2 1 

112 0 

12 2 1 

13 2 1 

2 

3 Z 

1 

3 

3 2 1 

5 

3 

2 1 

15.1 

5 3 2 

1 6 3 2] 

L 7 3 2 

1 8 3 2 1 

9 

3 2 

1 

10 

1 3 0 

11 

3 

0 

15.1 



to 


12 3 0 


0 

0-1.0 

0 

0-0.896 

0 

0 0.866 

0 

0 0.866 

0 

0-0. B66 

0 

0 0.866 

1 

0-0.6*»6 

1 

0-0.066 

1 

0 0.C66 

1 

0-0. 866 

1 

0 0.866 

0 

0- 0. 866 

0 

0-2.012 

0 

1- l.O 

0 

1-1.0 

0 

0-1.732 

0 

0-0.25 

0 

l-l.O 

0 

l-l.O 

0 

0 2,232 

0 

0 0.75 

0 

l-l.O 

0 

l-l.O 




BODY 

110 211 311 120 

2 3 1 3 3 0 

0 0 ~ 1.0 0 0 - 1.0 

1 O-l.O 0 0-0.5 

MU IT 


i 

! 


I 


0 0-0. B66 

0 0-0. B66 

1 O-l.O 

1 0 0.8b6 

I 0-0.866 

0 l-l.O 

0 0 0.5 

0 1 - 1,0 

0 0 1,732 

0 1 - 1.0 

0 0-1.232 


2 2 1 3 2 1 

I 0-1,0 

0 0 - 2.0 


15.1 

15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 
15.2 

15.0 

1 3 1 15.1 

15.1 

15.2 
15.2 


■ 



PDOGHAM L216A1 VERSION APR 2B, 77 NOW RUNNING. 
THE program is pari OF THE DVLOFL* SYSTEM 
DEVELOPED FOR NASA UNDER CONTRACT NASI-1391B. 
DATE Of RUN IS 77/09/17. 

TINE OF RUN IS U.1S.19. 


7.3 OUTPUT LISTING 



■*3 

4k. 


1 iru3L*I 

IlMUe PFOC^AN L216 sample PROfiLEM RING. STRUT, NACELLE, FUSELAGE 
( tUTLE 

(iriTLE TWELVE panels ( 2-R 1 NG , 1 - S TRU T ,6-NACE L L F INTERFERENCE, 
lilllLE 3-FUSELACE I N TE R F t RE NC E ) 

) ATITLE 

liTIILE three slender BODIES ( Y- NACE L LE , i-NAC E LL E , FUS EL ACE I 
( IT 1 TlE 

IATITlE seventeen strips (A-WlNG,l-SIRUT,6-NACeLLE INTERFERENCE, 

ia’itle 6-fusllace interference I 

I IM TL E 

tlTITLE THIRTY FOUR BOXES ( «- W I NG , 2- S TR UT , I i-N ACELLE INTERFERENCE, 
IATITlE 12-fuselage INTERFERENCE! 

t ATI U E 

IATITLE nine line elements (3-EACH SLENDER BODY! 

I AT ITlE 

(ATITlE three DOF IPLUNGE ,P1 TCM.ROLL/ 

I AI I TlE 

(ATITLE pitch ref X-0, ROLL REF V m Q, ^ > Q 
( All TlE 

iaiitle polynomial mode input 

« ATITLE 

(CASE 1 1 


CPU time I SECONDS! - .155 


INPUT parameters 


MACH number 
beta 

REFERENCE ChUBO 
REFERENCE SEMI-SPAN 

reference area 

REDUCED FFFOUENCIES 

total number of panels 12 

tutal nuhslp of bodies 3 

total number of strips it 

NUMBER OF VERTICAL PANEL STRIPS ON Y»0 PLANE 0 

NUMBER OF VERTICAL PANEL BOXES ON V-U PLANE 0 

motion is SYMMETRIC 


.8500 

.52t>8 

1.5000 

1.0000 

6.A000 

0 . 


.SOOOE^OO 


GEOMETRY DATA WILL BE SAVED ON TAPE SAVGEO 



AERO 


DATA hUL BE SAVED ON TAPE SAVAER 



■SJ 

A 


GEOMETRY ORTR 


PANEL GEOMETRY 


XI 


INBOARD 

X2 


1 nbuard 

X3 


OUTBOARD 

XA 


outboard 

VI 


inboard 

Y2 


OUTBOARD 

21 


INBOARD 

22 


OUTBOARD 


LEADING EDGE 
trailing EOCE 
LEADING EDGE 
TRAILING EDGE 

V CQDhOlNATE 

Y COORDINATE 

I coordinate 

I COORDINATE 


X COORDINATE 
X COORDINATE 
X COORDINATE 
X COORDINATE 



PANEL NO. 1 

INPUT 

VALUES 





XI 

- 2.000000 

X2 - 

A. 000000 

VI 

- .866000 

21 - 

.300000 

X3 

- 2.SOOOOO 

XA - 

A. 000000 

Y2 

« 2.000000 

22 * 

.500000 

3 

ChORONISE DIVISION 

BOUNDARIES 

FOR 

PANEL 1 IN FRACTION 

OF CHORD 


0. .5000EAO0 .1000E»0L 

3 SPANUISE DIVISION BOUNDARIES FOR PANEL 1 IN FRACTION OF SPAN 
0. .5000E*00 .1000E*01 


COORDINATES OF BOX CORNERS 

9 BOX CORNERS FOR PANEL I 


X 


Y 


I 



.20000E*01 

. S6600E«00 

.SOOQOEtOO 

.300C0fc«01 

. 866C0E»00 

.50000E»00 

.<>OOOOE«Ol 

.B6A00E»00 

.S00O0E»0O 

,22500E*0l 

. IA330E*01 

.SO0OOE*OO 

.312SOE*Ol 

.1A330E*01 

,S0O0OE»OO 

.AOOOOE+Ol 

.l<i330E»01 

.S0000£*00 

.25000E«01 

,20000E»01 

,50000E*00 

.32500E+01 

.200C0E«01 

.50000E*00 

.^OOOOE«01 

.20000Et0I 

.SOOOOEtOO 


STRIP L.E. CENTERLINE X-COORDINATES 
.2l250E*0l ,21750£»0l 

AVERAGE CHOPOLENGTH OF STRIPS 
.18750E*01 .162SOE*01 


—* PANEL NO. 2 INPUT VALUES ••• 


XI • 

2.500000 

X2 - 

A. 000000 

VI - 

2.000000 

21 - 

.500000 

X3 - 

3.000000 

XA - 

A. 000000 

V2 - 

3.000000 

22 • 

1.000000 


3 

CHOP ON I SE 

DIVISION BOUNOARI 

ES FOR PANEL 

2 

IN 

FRACTION 

OF 

CHORD 


0. 

.5C0CE»00 

.lOOOEAOl 






3 

SPANHISE 

DIVISION BOUNOARI 

ES FOR PANEL 

2 

IN 

FRACTION 

OF 

SPAN 


0. 

.SOOOEaOO 

.1000E*^0I 







COORDINATES OF BOX CORNERS 


9 BOX 

CORNERS FOR 

PANEL 2 

X 

V 

1 

. 25000E*0I 

•20000F*0l 

.50000E-00 

. 32500E*0l 

• 20UOO£*01 

•50000F»00 

. AO0(1OE*0l 

.200C0F*01 

.50000E*00 

. 27500E*0l 

.25000EA01 

.75000E*00 

. 33750£»0l 

.25000E*01 

. 7bOOOE+00 

.AOOCOE*01 

•25000e*01 

. 75000E»00 

.30000E*OI 

.30000E*01 

.10000E«-Ot 

.35000EA01 

.30000E»01 

.lOOOOE^Ol 




03 


.<.O0OO£ <01 . 300C0E<01 . 10000E*01 

STPIP L.E. CENTERLINE X-C OORD I H A T E S 
.262i0E»01 .28750C»01 

AVERAGE CHDBDLENGTH Of STRIPS 
. 13)S0E<0I . 112S0E<0I 


panel ho. 3 

INPUT 

VALUES •*» 





Xl - 2.S00000 

A2 - 

<1.000000 

VI - 

2.000000 

1 1 • 

.^00000 

X3 • 2.500000 

XV - 

N.OOOOOO 

V2 - 

2.000000 

22 - 

D.OOOOOO 


3 

ChDRDW I SE 

OIVISION RDUNOARJES FDR 

PANEL 

3 1 N 

FRACT [UN 

OF 

Chord 


0. 

.SOOOE*D0 , IOOOE*01 





2 

$PANn7S£ 

07VIS20N OOUNOARIfS FOR 

PANEL 

3 (N 

FPACT JON 

OF 

SPAN 


0. 

. 1000E<01 







COUhOlNATES OF BOA CORNERS 

6 eOX CORNERS FOR PANEL } 


X 

Y 

I 

250QDE <01 

. 20G0QE <01 

.50000E<00 

J2500E <01 

. iOOOUi‘01 

. jo.9ooe<oo 

SOCCOE <01 

. 20000E <01 

.50000E*00 

2SOOOE *01 

. 2000QF <01 

0. 

32SOOE <01 

•200C0E<01 

0. 

A00C0F<01 

.20000£*0l 

0. 


STRIP l.E. centerline X-COORO I NA TE S 
.250COE <01 

AVERAGE CMOROLEnCTh OF STRIPS 


lSOOOE*al 



••• PANEL NU. ^ INPUT VALUES 


XI ■ 

2.SOOOOO 

X2 ■ 

A. 000000 

ri - 

2.000000 

ll - 

0.000000 

X3 - 

z.sooooo 

■ 

X 

A. 000000 

V2 - 

2.A33000 

Z2 - 

-.2SOOOO 


3 CHOROWISE DIVISION BOUNDARIES FOR PANEL IN FRACTION OF CHORD 
0. .SOO0E»0O .lOOOEtOl 

2 SPANHISE DIVISION BOUNDARIES FOR PANEL A IN FRACTION OF SPAN 
0. .IOOOE*01 

COORDINATES OF BOX CORNERS 

b BOX CORNERS FOR PANEL A 


X 

Y 

2 

.2SOOOE«Ol 

.2OOCOE«0l 

0. 

. 32SOOE»01 

.20C00E»01 

0. 

. ACQ00E»01 

.20000E»01 

0. 

.2S000E*01 

.2A330E«01 

-.25000E*00 

.32SOOE*01 

.2A33OE»01 

-.ZSOOOE^OO 

.AOOOOE«01 

.2A330E*01 

-.25000E«00 


strip L.E. CENTERLINE X-COOROINATES 
.2SOOOE*OI 

AVERAGE CHOROlENCTH OF STRIPS 
. 1SOOOE«01 

PANEL NO. 5 INPUT VALUES ••• 


XI - 

2.SOOOOO 

X2 ■ 

A. 000000 

VI - 

2.A33000 

II - 

-.250000 

X3 - 

2.500000 

XA • 

A. 000000 

Y2 - 

2.A33000 

22 - 

-.T50000 


3 CHOROWISE DIVISION BOUNDARIES FOR PANEL S IN FRACTION OF CHORD 




S0C0E»OO 


1000E»0l 



2 SP*NMI$E DIVISION BOUNOADIES FOP PANEL S IN FRACTION OF SPAN 


0. .lOOOE^Ol 


COORDINATES OF BOX CORNERS 

6 BOX CORNERS FOR PANEL S 


X 

.Z50C0E*01 
. 325C0C*01 
.<tOOOOE«Ol 
.250006*01 
.325006*01 
.400006*01 


Y 

.243306*01 

.243306*01 

.243306*01 

.243306*01 

.243306*01 

.243306*01 


2 

-.250006*00 

-.250006*00 

-.250006*00 

75000E*00 

750006*00 

-.750006*00 


STRIP L.E. CENTERLINE X-COOROINATES 


. 250006 *01 

AVERAGE CHOHOLENGTH OF STRIPS 
. 150006*01 


•** 

PANEL NO. 

6 

INPUT 

VALUES ••• 

XI - 

2.500000 

X2 - 

4.0COOOO 

X3 - 

2.500000 

X4 - 

4.000000 

3 

CHOPDUlSE 

DIVISION 

BOUNDARIES 


0. 


.50006*00 

2 

SPANWISE 

DIVISION 

BOUNDARIES 


0. 


.10006*01 


VI - 

2.433000 

21 • 

-.750000 

V2 ■ 

2.000000 

22 - 

-1.000000 


FOR PANEL 6 IN FRACTION OF CHORD 
10006*01 

FOR PANEL 6 IN FRACTION OF SPAN 


COORDINATES OF BOX CORNERS 

6 BOX CORNERS FOR PANEL 6 



X 

V 

1 

.25000£*01 

.2<>330E« 01 

-. 75000E*00 

.32S00E*01 

.24330E*01 

-. 75000E*00 

.A0C0CE*01 

.2<>330E»0l 

-.75000E+00 

. 29000E*01 

.20000EAOI 

-.IOOOOE*OI 

.32SO0Et01 

,20000E»OI 

-.10000E«^01 

.bOOOOE»01 

.20000E»01 

-.10000E»01 


STRIP L.E. CEMTERLINE X-COOROINATES 
.2SOOOE»01 

AVERAGE CHOROLENGTH OF STRIPS 
. ISOOOE^Ol 



PANEL NO. 7 

INPUT 

VALUES ••• 

XI 

- 2.900000 

X2 ■ 

A. 000000 

X} 

- 2.900000 

XA - 

A.OOOOOO 


VI - 

2.000000 

Zl • 

-1.000000 

Y2 - 

1. 567000 

Z2 • 

-.750000 


i 

CHORDUISE 

DIVISION BOUNDARIES FOR 

PANEL 

7 IN 

FRACTION 

OF 

CHORD 


0. 

.SOOOE^OO .1000E*01 





2 

SPANN ISE 

DIVISION BOUNDARIES FOR 

PANEL 

7 IN 

fraction 

OF 

SPAN 


0. 

.lOOOE^Ol 







COOROIMATES OF BOX CORNERS 

b BOX CORNERS FOR PANEL 7 


X 

V 

Z 

. 25000E»01 

. 20000E+01 

-. 10000E*01 

.325rCE»01 

.2COOOE*Ol 

-. 100006*01 

. AOOOOE *01 

. 20000E*01 

-.100006*01 

.25C00E»01 

.15b70E>01 

-. 750006*00 

.32500E*01 

.15670F»01 

-.TSOOOEtOO 

.AOCOOEFOl 

.19b70E*01 

-.750006*00 


STRIP L.E. CENTERLINE K-COOROINATES 
.29000E«01 


00 



AVEKAGE CHOROLEM&TH OF STRIPS 


.ISOOOE^Ol 


••• PANEL NO. a INPUT VALUES 


Xl • 

2.500000 

X2 - 

4.000000 

VI - 

1.567000 

21 - 

-.750000 

X3 • 

2.500000 

XA - 

4.000000 

V2 • 

1.567000 

22 - 

-.250000 


3 CHOROUISE DIVISION BOUNDARIES FOR PANEL B IN FRACTION OF CHORD 
0. .SOOOEtOO .1000E«01 

2 SPANUISE DIVISION BOUNDARIES FOR PANEL B IN FRACTION OF SPAN 
0. .1000E»0l 


CUOFDINATES OF BOX CORNERS 

6 BOX CORNERS FOR PANEL 8 


X 

• 25OG0E >01 
.32S0OE»Ol 
.AOOOOEtOl 
.2SOOOE»01 
.32S00E*01 
.40COOE*01 


y 

.ISbTOE^Ol 
. 15670E*01 
. lSbT0E»01 
. 156T0E»0L 
.1S6T0E»01 
. tS670E»01 


I 

-.75000E+00 

750008*00 

75000E+00 

-.25000E+00 

-•25000E*00 

-.2SO0OE*OO 


STRIP L.E. CENTERLINE X-COORDINATES 
,25000E*01 


AVERAGE CHOROLENGTH OF STRIPS 
. ISOOOEvOl 


• •• 

PANEL NO. 9 

INPUT 

VALUES ••• 





XI 

• 2.500000 

X2 - 

4.000000 

Vi ■ 

1.567000 

21 - 

-.250000 

X3 

- 2.500000 

X4 - 

4.000000 

V2 «• 

2.000000 

22 • 

0.000000 



I 


3 CHCmOUISE DIVISION BOUNDAPIES FOR PRNEL 9 IN FRACTION OF CHORD 
0. .SOOOEtOO .10Q0E«01 

2 SPANUISE DIVISION BOUNDARIES FOR PANEL 9 IN FRACTION OF SPAN 
0. .lOOOE^Ol 


COORDINATES OF BOX CORNERS 

6 BOX CORNERS FOR PANEL 9 


X 

,25000E»01 
. 32500E*01 
.AOOOOE*Ol 
.2S000E »01 
.32SOOE«Ol 
. AOOOOE*01 


V 

.15670E»01 
. 1S670E*01 
.15670001 
.20000€*0l 
.2C000E+D1 
.20000E*01 


z 

-.25000E*00 

-.25000E+00 

-.25000E«00 

0 . 

0 . 

0 . 


STRIP L.E. Centerline x-coordinate$ 
.25000E«01 


AVERAGE CHORDLENGTH OF STRIPS 
. ISOOOEAOl 


PANEL NO. 10 INPUT VALUES 


XI - 

2.000000 

X2 - 

A.OCOOOO 

Yl • 

0.000000 

Zl ■ 

1.000000 

X3 - 

2.000000 

X6 ■ 

A . 000000 

V2 • 

•fi6fr000 

22 • 

.500000 


3 

CHOROUISE 

DIVISION BOUNDARIES 

FOR 

PANEL 

10 

IN 

FRACTION 

OF 

CHORD 


0. 

.50COEAOO 

1000E»01 






3 

SPANN I SE 

DIVISION BOUNDARIES 

FOR 

panel 

10 

IN 

fract ion 

OF 

SPAN 


0 . 


SOCOE^OO 


lOOOEROl 



CODPOINATES Of BOX CODNERS 


9 BOX CORNERS FOR PANEL 10 


X 

.20000E »01 
. 3COOOF *01 
.LOOOOe »01 
. ZOCOOE^Ol 

. 30coce*oi 

• AOOCOE *01 
.20000E*0l 
.30000Et01 
. AOOOOE^Ol 


Y 

0 . 

0 . 

0 . 

.A33COE»00 
.A330OE^00 
.A3300E*00 
. 86600E*00 
.UbAOOEtOO 
.8&600E400 


I 

.lOOOOE»Ol 
.lOOOOE^Ol 
.10000E*01 
,75000E*00 
. 75000£*00 
. 7S00OE+00 
.SOOOOE^-OO 
.SOODOEi-OO 
.SOO0OE*OO 


STRIP L.E. CENTERLINE X-COOROINATES 
.2C0O0E«Ol .ZOOOOE^Ol 

AVERAGE CHOROLENGTH OF STRIPS 
.20000E<'01 .20000EF01 


panel no. 11 INPUT VALUES **• 


XI - 

2.000000 

X2 - 

4». 000000 

Yl - 

.866000 

21 • 

.sooooo 

X3 - 

2.000000 

» 

A. 000000 

Y2 ■ 

.866000 

22 • 

-.sooooo 


3 CHOROrilSE DIVISION BOUNDARIES FOR PANEL U IN FRACTION OF CHORD 
0. .SOOOE^OO .iOOOExOl 

3 SPANMlSE DIVISION BOUNDARIES FOR PANEL 11 IN FRACTION OF SPAN 
0. .S000E«^00 .LOOOE^Ot 


COORDINATES OF BOX CORNERS 

9 BOX CORNERS FOR PANEL 11 


V 


I 



.200UOE«01 
,3COCOb*OI 
. ^OOOOEtOl 
.20000F»01 
.30000E*01 
,40000F»01 
.2000UE*0t 
.30000E*01 
.^OOOOE«Ot 


. H6b0UE«00 
.B66COE*00 
. 86o00E»00 
.BC.bOOE»00 
.866C0E*00 
.86b00F*00 
. 8b600E»00 
.86&r0E»00 
.BbbOOEtOO 


.50Q00F«00 

.50000E*0a 

.50000E*00 

0 . 

0 . 

0 . 

-.SOOOOE«00 

-.50000E*00 

-.30000£»00 


SThlP L.E. CENTERLINE N-COOROINATES 
.20000E«0L .ZOOOOE«01 

AVERAGE CHOROLENGTH OF STRIPS 
.20000E*01 .20000E«01 

••• PANEL NU. 12 INPUT VALUES 


XI > 

2.000000 

X2 - 

4.000000 

VI - 

.366000 21 ■ 

-.500000 

X3 - 

2.000000 

X4 - 

4.000000 

V2 • 

0.000000 22 • 

-1.000000 


3 CHOROWISE DIVISION BOUNDARIES FOR PANEL 12 IN FRACTION OF CHORD 
0. .S000E»00 .lOOOE+Ol 

3 SPANMISE DIVISION BOUNDARIES FOR PANEL 12 IN FRACTION OF SPAN 
0. .5000E»00 .lOOOEAOl 


COORDINATES OF BOX CORNERS 

9 BOX CORNERS FOR PANEL 12 


X 

.20000E«01 
.300COE*01 
. A00nOE«Ol 
. 20000E*01 
. 30000E*01 
•AOOOOEtOl 
.20000EAOI 
.30000E*01 
,AUOOOE*Ot 


V 

. BbbOUEFOO 
.Rb600E*00 
.SbbOOE^OO 
.A33C0E*00 
.43300E*00 
.43300E«00 
0 . 

0 . 

0 . 


2 

-.5')000E»00 

~.S000OE»00 

-.500006*00 

-.T5000e*00 

-.75000E*00 

-.750006*00 

100006*01 

-.10000E*01 

-.100006*01 



STRIP L.E. CENTCRLINE X-CUURO J NA TES 
.20000E»01 .20000E*01 

AVERAGE CHaROLENGTH OF STRIPS 


ZOOOOEFOl 


20000E»01 



• PANEL OAtA • 


ARKAV OF CUMULATIVE tlUX NUMBERS FOR ALL PANELS 

4 8 10 12 1<. lb 18 20 22 2b 

DIHEDRAL ancles FOR ALL PANELS 

0.00 2b.S7 270.00 330,00 270.00 

150.00 90.00 30.00 330.00 270.00 



30 3A 


210.00 

210.00 



00 

00 


• STRIP DATA • 


CENTERLINF 


V 


z 

OELVS 

OELZS 

. 1 UR50E 

*01 

•SOOOOOE*00 

-567000E*00 

0. 

. 1 71LS0F 

♦ 01 

.500000E*00 

.5670C0E*00 

0. 

.22SCLCE 

♦ 01 

• 625000E *00 

.5000COE*00 

.250000E*00 

.27SUUDE 

♦01 

• a 75POOE *00 

.5000001 *00 

,25nD00E*00 

.200C00F 

♦ 01 

.2!)001:CE*00 

0. 

.50000CE*00 

.221 OtOF 

*01 

-.12S00HE*00 

.A33000t*00 

.2500009*00 

.2A1 jOOF 

♦01 

-.sooooor+co 

0. 

• somnoF *00 

.2216^0*^ 

*01 

-.B75CP0E*00 

.A33000F*00 

.250000E+00 

. 17B3SQF 

♦ 01 

-.87SO0OE*O0 

.A33000E*00 

.2500009*00 

.156T00E 

♦ 01 

-.SOOOCOE*00 

0. 

.500000E*00 

. 1783SOE 

♦ 01 

1250COE*00 

.A33000F*00 

.2500UOE*00 

•216S00E 

♦ 00 

.a75000F*00 

.a33000E*00 

.2500009*00 

.6ANS00C 

♦ 00 

• 625000E *C0 

. A33000E+00 

.2500009*00 

.86<>CCQE 

♦ 00 

.2SQ0Q0E*00 

Q. 

.5000001*00 

. SbbOOOE 

♦ 00 

-.2500009*00 

0. 

.500000E*00 

.6A9500E 

+ 00 

-.625000E*00 

.A33000E+00 

. 250000E*00 

.216500E 

♦ 00 

-.87SO00E*00 

.A33000E*00 

. 250000E*00 


I 

I 


fpact ion of 
PANEL SPAN 

.250C00E *00 
. TSOOOOEtOO 
.2S0000F *00 
. T50000E *U0 
.SOOOOOEtOO 
.bU000OE*0O 
.500000E »00 
.500000F «00 
.5DOOOOE *00 
.500000E *00 
,SQ0OO0E*O0 
.ZSOOOOf *00 
. 75000un*00 
•2S0OOUE»OO 
.TSOOOOEtOO 
.250000£*00 
.TS0000E*00 


• BOX GEOMETRY* 


PART I 

COOROINATES OF BOX EDGES AT 1/A CHORD POINTS 


INBOARD 


XVI 


.22S0E*01 

.SC.60E »00 

. snooE»oo 

.12S0E*01 

.6b60E »0O 

.5000E»00 

.2A6SE*01 

. IR33E»01 

.5000£»00 

.33AAE*01 

. IA33F»01 

.5000E*00 

.26BaE*Ot 

.2O0OE»0l 

. 5000P*00 

. 3A36E*01 

.2000E«01 

. S0O0E»O0 

.2906E»D1 

•2500F*01 

. 7S0OE»OO 

.3S31E*Ol 

.250UE lOl 

. 75C0E*00 

.26HB£»01 

,2000f »01 

.5000E»00 

. 3A3BE«Ol 

.20C0£*0l 

.50008*00 

.266BE»ni 

•20COF»Ol 

0. 

. 3A3BE*01 

.2000E»01 

0. 

.2ASBE*Ot 

.2A33F*01 

-.2300E*00 

.3A3BE»01 

.2A33E*01 

-.2500E *C0 

.26SaE*01 

.2A33E*01 

-.7500fc*00 

.3A3BE*Ol 

.2r33E«01 

-.7500F»00 

.2»8BE*01 

.2000E*01 

1000E*0i 

.3A3BE*01 

.2OCOE»01 

-.10008*01 

.26B8E*01 

.1567F»01 

7500t*00 

.3A3B£*01 

.136TE»01 

7500t*00 

.268B£»01 

.15fcTF»Ol 

-.2500E*C0 

.3A3e£*01 

. 1567E*01 

-.25007*00 

-225CE*0l 

0. 

. lC00k*01 

.32S0E*Ol 

0. 

.10008*01 

.2230E*01 

.A330E»00 

.75008*00 

.32S0b*0t 

.A33CE*00 

. 7500E*C0 

.22SOE*Ol 

.8660E*00 

.50008*00 

.32»0E*01 

.SbbOE^OO 

.5000E*00 

.2250E*0I 

.fl660E«G0 

0. 

.32S0E*01 

.6b60E*00 

0. 

.22S0E*01 

.B6b0E*00 

-.50008*00 

.3250E*01 

.Bbb0E«^00 

-.50008*00 

•22S0E*01 

.A330E*00 

-.7500E*C0 

.32S0E*Ol 

.A330E*00 

-.7500E*00 



OUTBOARD 


BOX 

X 

Y 

1 

NO. 

.2AbSE*01 

. 1633E*01 

.50006*00 

1 

.33S<,F*01 

.16338*01 

.50008*00 

2 

.26808*01 

.20006*01 

.50006*00 

3 

.3R3BE*01 

.20008*01 

.50006*00 

A 

.290bE*0l 

.250CE*01 

.75006*00 

5 

.3531E*0l 

.25006*01 

.75006*00 

6 

.3125E»0t 

.30006*01 

.10006*01 

7 

.3b25E»01 

.30006*01 

.10006*01 

8 

. 2bb8E*01 

.20006*01 

0. 

9 

.3A3BE*01 

.2000E*01 

0. 

10 

.26B8E*01 

.2633E*01 

-.25006*00 

11 

.36388*01 

.26336*01 

-.25008*00 

12 

.2bB8E*0l 

.26336*01 

-.75006*00 

13 

.3638E*01 

.26336*01 

-.75006*00 

lA 

.2680E*01 

,20006*01 

-. 1000E»01 

15 

.363flF*0l 

.20006*01 

-.10006*01 

16 

.2688E*01 

.15676*01 

-.75006*00 

17 

.36388*01 

.15678*01 

-.75008*00 

18 

.26BBE*01 

. 15676*01 

-.25006*00 

19 

.36388*01 

. 15676*01 

-.25006*00 

20 

.26888*01 

.2000F *01 

0. 

21 

.3638E*01 

.20006*01 

0. 

22 

.22508*01 

.63308*00 

.75006*00 

23 

.32508*01 

.63306*00 

.75006*00 

2A 

,2250F*0l 

.86606*00 

.50006*00 

25 

.32508*01 

,86606*00 

.50006*00 

26 

.22508*01 

. 86608*00 

0. 

27 

.3250E*01 

. 86606*00 

0. 

28 

.22506*01 

.66606*00 

-.50006*00 

29 

.3250E*01 

.86606*00 

-.50006*00 

30 

.22508*01 

.63306*00 

75006*00 

31 

.32 50E*01 

.63306*00 

-.75006*00 

32 

.22508*01 

0. 

-.10006*01 

33 

.3250E*0t 

0. 

-.10006*01 

3A 


8 


PART 2 

COORDINATES OF RUX SENDING AND RECEIVING POINTS 


1/4 

CHORD POINTS 

3/9 

CHORD POINTS 


BOX 

X 

V 

1 

X 

Y 

1 

NO. 

.235<>E»01 

.11506*01 

.500r£*00 

.28286*01 

, 1 1506*01 

.50006*00 

1 

•J297E»Ol 

.11506*01 

.50006*00 

.37666*01 

,11506*01 

.50006*00 

2 

,2578E*0I 

. I717E*01 

.5000E*00 

.29896*01 

.17176*01 

.50006*00 

3 

.3391£«Ol 

. 1717E*01 

.500CE*00 

. 37976*01 

.17176*01 

.50006*00 

9 

.2797E*01 

.22506*01 

•6250E*00 

.31916*01 

. 22506*01 

.62506*00 

5 

,3‘i84E*0l 

.2250E*01 

.62506*00 

. 3H28E *01 

.22506*01 

.62506*00 

6 

. 301ibE>01 

.2750E*01 

•87506*C0 

.32976*01 

.27506*01 

.87506*00 

T 

. i5T8E*0V 

.27 506*01 

.97506*00 

.38596*01 

. 27506*01 

.87506*00 

8 

,2688E*C1 

.20006*01 

•2500£*00 

,30636*01 

.20006*01 

.25006*00 

9 

. 3438c»01 

, 200CI *Cl 

.25006*00 

.38136*01 

.20006*01 

.25006*00 

10 

.268BE*C1 

.2217E*01 

12506*00 

.30636*01 

.22176*01 

-.12506*00 

11 

. 3438E«01 

.22176*01 

-.12506*60 

.38136*01 

.22176*01 

-.12506*00 

12 

.2680E»OI 

. 29336*01 

-.50006*00 

. 30636*01 

.29336*01 

-.50006*00 

13 

.343BE»01 

.29336*01 

-.50006*00 

. 38136*01 

.29336*01 

-.50006*00 

19 

.2688E>0t 

.22176*01 

-.8750£*00 

. 30636*01 

.22176*01 

-.87506*00 

IS 

,343aE*Ol 

.22176*01 

-.8750EtOO 

.38136*01 

.22176*01 

-.87506*00 

16 

. 2t>88E*01 

.17896*01 

-.8750E *C0 

. 30636*01 

.17896*01 

-.87506*00 

1 7 

. 3438E«CI 

.17896*01 

-.87506*00 

.38136*01 

.17896*01 

-.87506*00 

18 

.268BE »01 

. 15676*01 

50006 *00 

.30636*01 

.15676*01 

-.50006*00 

19 

. 3<r38Et01 

, 1567E*0I 

-. 5000E *00 

,38136*01 

. 15676*01 

-.50006*00 

20 

. 2b38E»01 

.17896*01 

-.12506*00 

.30636*01 

.17896*01 

-.12506*00 

21 

.3438t*0l 

. I789£*0l 

-.12506*00 

.38136*01 

.17896*01 

-.12506*00 

22 

.22506*01 

.21t5E*rO 

. S750E *00 

.27506+01 

.21656*00 

,87506*00 

23 

. 3t50£*01 

.2U5F*00 

.87506*00 

.37506*01 

.21656*00 

.87506*00 

29 

. 22506*01 

.69956*00 

.62506*00 

,27506*01 

.69956*00 

.62506*00 

25 

.32506*01 

.6'.996*C0 

.62506* no 

.3760E+01 

.69956*00 

.62506*00 

26 

.2250£*01 

.8660E*00 

.2500E*00 

.27506*01 

. 86606*00 

.25006*00 

27 

.32506*01 

. A660E*00 

.25006*00 

.37506*01 

. 86606*00 

.25006*00 

28 

.22506*01 

.86oC6*00 

-.25006*00 

,27506 *01 

. 86606*00 

-.25006*00 

29 

.32506*01 

.86606*00 

-. 2500E*00 

.37506*01 

. 86606*00 

-.25006*00 

30 

.22506*01 

.6945E*00 

-.62506*00 

.27506*01 

.69956*00 

-.62506*00 

31 

.32506*01 

.69956*00 

-.62506*00 

.37506*01 

.69956*00 

-.62506*00 

32 

.22506*01 

.21656*00 

-.87506*00 

.27506*01 

.21656*00 

-.87506*00 

33 

.32506*01 

.21656*00 

-.87506*00 

.37506*01 

.21656*00 

-.87506*00 

39 



X-COOBOINATES OF SENDING POINTS AS FRACTIONS OF STRIP CHOROLENCTH 


. 125000E *00 
.62SCl'0E *C0 
.I2SCO0E*OO 
.62SnC0E*C0 
. 125000E*00 
.62SO00E*00 
. 12SOUOE *00 


.625000F*00 

.125000E*00 

.62S003F*00 

.12S0C0£*00 

.6250006*00 

.1250006*00 

.6250006*00 


.1250006*00 
.6250006*00 
.1250006*00 
.6250C0E*00 
. 1250006*00 
.6250006*00 
.1250006*00 


.6250006*00 

.125001)6*00 

.6250006*00 

.1250006*00 

.6250006*00 

.1250006*00 

.6250006*00 


CENTERLINE ChORDLENGTH UF BOXES 


.937500E*00 
.6(175006 *00 
. 75OCO0E *00 
. 750000F *00 
.7500006*00 
.1000006*01 
.1000006*01 


.9375006*00 
.5 6 2 5 0 06 * 00 
.7500006*00 
. 7500006*00 
. 7500006*00 
.1000006*01 
.1000006*01 


.8125006*00 
.5625006*00 
. 750000F*00 
. 7500006*00 
. 100000E*01 
.1000006*01 
.1000006*01 


.ai2500£*00 

.7500006*00 

.7500006*00 

.7500006*00 

.1000006*01 

.1000006*01 

.1000006*01 


BOX widths 


.5670006 *00 
.55901 7E*00 
.5999H9E *00 
.6999896 *00 
.6999896*00 
.6999896*00 
.699989E*00 


.567000F*00 
.5590176*00 
.■<.9 99 89^*00 
.9999896*00 
.9999896*00 
.5000006*00 
.9999896*00 


.5670006*00 

.5590176*00 

.5000006*00 

.9999896*00 

.9999896*00 

.5000006*00 

.9999896*00 


.5670006+00 

.5000006*00 

.5000006*00 

.5000006*00 

.9999896*00 

.5000006*00 

.9999896*00 


.1250006*00 

.6250006*00 

.1250006*00 

.6250006*00 

.1250006*00 

.6250006*00 


.6875006*00 
.7500006*00 
.7500006*00 
.7500006*00 
.1000006*01 
. 1000006*01 


.5590176*00 

.5000006*00 

.9999896*00 

.5000006*00 

.9999896*00 

.5000006*00 


BUOY geometry 


••• BODY NO. 1 INPUT VALUES 

CENTER OF BODY CDOROI NATES Y - 2.000000 2 • -.SOOOOO 

YFLAC > 1 2FLAG - 0 MODE SHAPE COEFFICIENT • 1.000000 

BODY BOA LIMITS ON INTERFERENCE PANELS FOR THIS BODY a 11 22 

A body element endpoints for body no. I 

200000EF01 .2SOOOOE«01 .32S000E+01 .ASOOOOE*Ol 

A BODY RADII AT ELEMENT ENDPOINTS FOR BODY NO. 1 
-SOQOCOEaOO .SOOOCOEfOO 0. 

••• BODY NO. 2 input VALUES 

CENTER OF BUOY COORDINATES Y « 2.000000 I ■ -.500000 

YFLAG • 0 2FLAG - I MODE SHAPE COEFFICIENT • 1.000000 

BODY BOX LIMITS ON INTERFERENCE PANELS FOR THIS BODY > 11 22 

A BODY ELEMENT ENDPOINTS FOR BODY NO. 2 
200000E*’01 .250000E»01 .325000E»01 .ASOOOOEfOI 

A BODY RADII AT ELEMENT ENDPOINTS FUR BODY NO. 2 
.5000D0E*00 .500000E»00 0. 

••• BODY NO. 3 INPUT VALUES 

CENTER OF BODY COUROINATES Y > 0.000000 / • 0.000000 

YFLAC a 0 ZFLAC a l RODE SHAPE COEFFICIENT a 1.000000 

BODY BOX LIMITS ON INTERFERENCE PANELS FOR THIS BODY a 23 3A 

A BODY ELEMENT ENDPOINTS FOR BODY NO. 3 




0. .200000t*01 .'i000C0f»01 .600000E*01 

BODY B*DII AT ELEMEMI ENOPOJNTS FOR BODY NO. 3 
0. .100000E*01 .100000£*0l 0. 


CO 

03 


• BODY ABHAYS • 


receiving point 


SENDING 

LENGTH 

AVG. RADIUS 

80X 

X 

Y 

1 




N 

.22SOE»Ol 

.200QE*01 

- . 50006*00 

.22506*01 

. 50006*00 

.25006*00 

35 

.2B75E+01 

.20006*01 

-.50006*00 

.28756*01 

.75006*00 

.50006 *00 

36 

. 3875E*01 

,2000E*0l 

-.50006*00 

.38756*01 

.12506*01 

.25006*00 

37 

•2250EA01 

.20006*01 

-.50006*00 

.2250E*01 

. 50006*00 

.25006*00 

38 

.2875E*OX 

.2000E*01 

-.50006*00 

.23756*01 

.75006*00 

.50006*00 

39 

. 38756*01 

.20006*01 

-.50006*00 

. 38756*01 

.12506*01 

.25006*00 

*tO 

. 1000E*01 ' 

0. 

0. 

.10006*01 

.20006*01 

.50006*00 

Al 

.3000E*01 

0. 

0. 

.30006*01 

.20006*01 

.10006*01 

A2 

.5000E*01 

0. 

0. 

.50006*01 

.2000E*0l 

.50006*00 

A3 


END OF GEOMETRY DATA 


HOOAL DATA 


POLYNOHtAL DATA FOP MODES OEFINEO 

PPINI PPESSURES AND GEN. FORCES 
total number of MOOES 


total number of modes 3 

total number of coefficients AO 

NUMBER OF COEFFICIENTS FOR PANELS 3A 

NUMBER OF COEFFICIENTS FOR BODIES 6 


POLYNOMIAL MODE DEFINITION 


PANEL MODE POMER OF POWER OF 


COEFFICIENT 


NO. 


NO. (X / SI <Y / SI 


1 1 

2 1 

A 1 

6 1 

7 1 

5 1 

10 1 

12 1 

1 2 

2 2 

A 2 

6 2 

T 2 

9 2 

10 2 

12 2 

1 3 

1 3 

2 3 

2 3 

3 3 

3 3 

A 3 

A 3 

9 3 

5 3 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I 

0 

1 

0 

1 

0 

I 

0 

1 


- 1.0000 
- .89A0 
-•BbbO 
.ObbO 
. RbbO 
-.B6b0 
flobO 
.BbbO 
-I .0000 
-.09a0 
“.BbbO 
.BbbO 
. BbbO 
-.BbbO 
BbbO 
• BbbO 
-.BbbO 
- 1.0000 
- 2.0120 
-l.OOOO 
.5000 
- 1.0000 
-1.7320 
- 1.0000 
-.2500 
• 1.0000 


6 3 0 

6 3 0 

7 3 0 

7 3 0 

6 3 0 

8 3 0 

9 3 0 

9 3 0 

10 3 0 

10 3 0 

10 3 0 

10 3 0 

10 3 0 

10 3 0 

10 3 0 

ID 3 0 

10 3 0 

10 3 0 

10 3 0 


0 

1 

0 

1 

0 

1 

0 

I 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1.7370 
- 1.0000 
2.2320 
- 1.0000 
. 7500 
-l.OOOO 
-1.2320 
- 1,0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0030 
0.0000 
0.0000 
0.0000 
0.0000 


PQLVNOHIAL MODE DEFINITION 


PANEL MODE POwEA OF POltER OF COEFFICIENT 


NO. NO. (X /SI (Y / SI 


10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 

10 3 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

o 

0 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


the 39 8-MATRIX ELEMENTS FOR MODE NO. 1 


531563E»00 -,531563E»00 -.96068B£»00 -.960688E»00 -.393586E*O0 
3297I.3E»00 3297 93E ♦00 0. 0. .329793E^OO 
329793E»00 -.329793E»00 -.932990£^00 -.932990£^00 -.932990E+00 
432990E«00 .932990E^OO .932990E»00 .932990E^00 


-. 3935866+00 
.3297936+00 
-.932990E+00 


-.281 116E+00 
.329793E+00 

0 . 


-.281116E + 00 0, 
.329T93E+00 0. 

0 . 0 . 


THE 39 B-MATRIX ELEMENTS FOR NODE NO 


2 



125416E*01 

-.0727<.6E*OO 

-.B727^6f*00 

.97<i229E*O0 


175250E«0l 
IIW 30E*01 
-.lllb30E»01 
.U0722E»01 


-.118771E*01 

0 . 

- .97't229E*00 
• 979229E»00 


15t.202£*01 

0 . 

-,U0722E«0l 
. l'tO722E*0l 


-.9(>09b 7E»00 
.0 72 7<,SE«-OO 
-.979229Et00 


-.119718E*01 
. 1 1 lb lOEt-01 
-.140722E*01 


-.697739E*00 1005876*01 0. 

.e72796E*00 .1115306*01 0. 

0 . 0 . 0 . 


THE 


39 0-MATRlX 


ELEMENTS FOB MODE NO. 


3 


-.bll031E*00 

-.7432326*00 

-.SS5736E*00 

0 . 


-.611031E*00 

-.743232£*00 

-,55573bE*00 

0 . 


-.7907706*00 
- . I075OOE*OO 

0 . 

0 . 


790770E*00 

-.1875006*00 

0 . 

0 . 


-.B80b82E*00 

.5557406*00 

0 . 


-.880682E*00 

.5557406*00 

0 . 


-.89b339E*00 

•74323bE*00 

0 . 


-.89b339E*00 

.7432386*00 

0 . 


,9375006-01 

.1875006*00 

0 . 


.93 7500 6-01 
.1875006*00 

0 . 


CO 

-7 


PULYNOHUL HOOE OEFlNniDN 


SODV 

MODE 

POWER 

OF 

POWER 

OF 

COEFFICIENT 

NO. 

NO. 

IX / 

Si 

(V / 

SI 

(ARQI 

2 

1 


0 


0 

-1.0000 

3 

1 


0 


0 

-1.0000 

2 

2 


1 


0 

-1.0000 

3 

2 


t 


0 

-1.0000 

1 

3 


0 


0 

-.5000 

2 

3 


0 


0 

-2.0000 


>.62$OOOEtOO -.lOOOOOE^Ol -,200000E»01 -.100000E»01 


THE 9 B-HATRtX ELEMENTS FOR HOOE NO. 

0. -.2S0000E*00 •'.7SOOOOE»O0 


the 

0 . 

THE 


9 B-MATR[X elements FOR MODE NO. 2 

-.Sb2S00E*00 -.21S62SE«^01 -.2A218BE»01 -.lOOOOOEfOl >.600000E«01 •>. SOOOOOEfOI 
9 B'MATRIX ELEMENTS FOR HOOE NO. 3 


12SOOOE*00 >.37S000E*00 -.312SOOEtOO -,S00000E»00 -.ISOOOOEtOl 


12SOOOE«01 0 


0. 


0 , 



• REOUCEO FREQUENCY • 0. 



THE 

3A00HNHA SH 

W 

ELEMENTS FOR MODE 

NO. 

1 



0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 


the 

3AD0WNUA SH 

w 

ELEMENTS FOR MODE 

NO. 

2 



.100000E»01 

0. 

.tOOOOOE»01 0. 


lOOOOOEvO L 

0. 

.89A0OOE»O0 

0. 

.894000E*00 0. 


89Ji000E»aQ 

0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 


THE 

3900UNUASH 

N 

ELEMENTS FOR MODE 

NO. 

3 



0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 



THE 

9 

HJ MATRIX ELEMENTS FOR 

MODE 

NO. 

1 

0. 

0. 

0. 


0, 

0. 



0. 

0. 

0, 

0. 


0, 

0* 



0. 



THE 

9 

HJ-PRIHE ELEMENTS FOR 

MODE 

NO. 

1 

0. 

0. 

0. 


0. 

0. 



0. 

0. 

0. 

0. 


0. 

0. 



0. 

COL. NO. L 

UF line 

ELEMENT PRESSURE COEFFICIENT 

MATRIX 



0. 

0. 

0. 


0. 

0. 



0. 

•0. 

0. 

0. 


0. 

0. 



0. 


The 9 UJ MATRIX ELEMENTS FOR MODE NQ. 2 


0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 


.lOOOOOE«Ol 0. .B9A000E*00 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 0 . 0 

0 . 0 . 




100 


0. 

0. 

0. 

0. 


0. 


0. 

. lOOOOOE^Ol 

0. 

.100000E»01 

0. 

.lOOOOOE^Ol 0. 

.lOOOOOEtOl 0. 


•lOOOOOEtOl 

0. 

.lOOOOOEtOl 

0. 





THE 9 

MJ-PPI ME 

ELEMENTS 

FOR MODE 

NO. 

2 





0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 



COL. NO. 

2 OF LINE 

ELEMENT PftESSOdE COEFFICIENT MATRIX 







0. 

0. 

0. 

0. 


0. 


0. 

•628319E*01 

0. 

0. 

0. 

-.251327E«01 0. 

.31LIS9E^01 0. 


0. 


0. 

-.3U159E*01 

0. 





THE 9 

WJ MATRIX 

ELEMENTS 

FOR MODE 

NO. 

3 





0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 





THE 9 

HJ-PRIME 

ELEMENTS 

FOR MODE 

NO. 

3 





0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 



COL. NO. 

3 OF line 

ELEMENT PRESSURE COEFFICIENT MATRIX 







0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 


0. 

0. 

0. 





HODE NO 1 


BOX PI^ESSURE DIFFERENCES 


NO. 

XOC 

X 

Y 

1 

PRESSURE 

REAL 

DIFFERENCE 

IHACINARV 

1 

.12500 

2. J5938 

1. 19950 

.50000 

0. 

0. 

2 

.62500 

3.29688 

1.19950 

.50000 

0. 

0. 

3 

.12500 

2.57813 

1.71650 

.50000 

0. 

0. 


.6 25 00 

3.39063 

1. 71650 

. 5000Q 

0. 

0. 

5 

. 1 2500 

2. 19688 

2.25000 

.62500 

0. 

0. 

6 

.62500 

3.98933 

2.25000 

.62500 

0. 

0. 

7 

.12500 

3.01563 

2.75000 

.87500 

0. 

0. 

B 

.6 2500 

3.57813 

2. 75000 

.87500 

0. 

0. 

9 

. 1 2500 

2.68750 

2.00000 

,25000 

0. 

0. 

10 

.62500 

3.93750 

2.00000 

.25000 

0. 

0. 

U 

.12500 

2.68750 ■ 

2.21650 

-.12500 

0. 

0. 

12 

,62500 

3.93750 

2.21650 

-.12500 

0. 

0. 

13 

.12500 

2.68750 

2.93300 

-.50000 

0. 

0. 

19 

.62500 

3.93750 

2.93300 

-.50000 

0. 

0. 

15 

.12500 

2. 68750 

2.21650 

-.87500 

0. 

0. 

16 

.6 25 00 

3.93750 

2.21650 

-.67500 

0. 

0. 

17 

.12500 

2.68750 

1. 78 J50 

-.87500 

0. 

0 . 

IB 

.62500 

3.93 750 

1 . 78350 

-.87500 

0. 

0. 

19 

.125C0 

2.68150 

1.56700 

-.50000 

0. 

0. 

20 

.62500 

3.93750 

1.56700 

-.50000 

0. 

0. 

21 

.12500 

2.68750 

1.78350 

-.12500 

0. 

0. 

22 

.62500 

3.93750 

1.78350 

-.12500 

0. 

0. 

23 

.12500 

2.25000 

.21650 

.87500 

0. 

0. ' 

29 

.62500 

3.25000 

.21650 

.87500 

0. 

0. 

25 

.12500 

2.25000 

.69950 

.6 2500 

0. 

0. 

26 

.62500 

3.25000 

.69950 

.62500 

0. 

0. 

27 

.12500 

2.25000 

. 86600 

.25000 

0. 

0. 

2U 

.62500 

3.25000 

. 86600 

.25000 

0. 

0. 

29 

.12500 

2.25000 

.86600 

-.25000 

0. 

0. 

30 

.62500 

3.25000' 

.86600 

-.25000 

0. 

0. 

31 

.12500 

2.25000 

.69950 

-.62500 

0. 

0. 

32 

.62500 

3.25000 

'.69950 

-.62500 

0. 

0. 

33 

.12500 

2.25000 

.21650 

-.87500 

0. 

0. 

39 

.62500 

3.25000 

.21650 

-.87500 

0. 

0. 


O 
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PANEL 

sectional OERIVATIVES 


NODE NO I 


STRIP NO. 

V 

1 

VOS 

LIFT COEFFICIENT 

MOMENT COEFFICIENT 





REAL 

IMAG 

REAL 

IMAG 

1 

1.1A9S 

.5 000 

1.1695 

0.000000 

O.OOQQOO 

0.000000 

0.000000 

2 

1.7165 

.5000 

1.7165 

0.000000 

0.000000 

0.000000 

0.000000 

i 

2.2500 

.6250 

2.2500 

0.000000 

0.000000 

o.oooroo 

0.000000 

A 

2.7500 

.8750 

2. 7500 

0.000000 

o.onoooo 

o.ooocoo 

0.000000 

S 

2.0000 

.2500 

2.0000 

0.000000 

0.000000 

O.OOOOOO 

0.000000 

6 

2.2165 

-.1250 

2.2165 

0.000000 

0.000000 

o.ooooco 

0.000000 

7 

2.A330 

-.5000 

2.6330 

0.000000 

0.000000 

0.000000 

0.000000 

8 

2-2165 

-.8750 

2.2165 

0.000000 

0.000000 

O.OOOOOO 

0.000000 

9 

1.7B35 

-.8750 

1.7835 

0.000000 

o.oocooo 

0.000000 

0.000000 

10 

1.5670 

-.5000 

1.5670 

0.000000 

0.000000 

0 . 000000 

0.000000 

11 

1.TH35 

-. 1250 

I.7B15 

0.000000 

0.000000 

o.ooooco 

0.000000 

12 

.2165 

. 6750 

.2165 

0.000000 

0.000000 

0.000000 

0.000000 

13 

.6A95 

.6250 

.6695 

0.000000 

0.000000 

0.000000 

0.000000 

lA 

.8660 

.2500 

.8660 

0.000000 

0.000000 

0.000000 

0.000000 

1$ 

.8660 

-.2500 

.8660 

0.000000 

0,000000 

0.000000 

0.000000 

16 

.6595 

-.6250 

.6695 

0.000000 

0.000000 

0.000000 

0.000000 

17 

.2165 

-.8750 

.2165 

0.000000 

0.000000 

0.000000 

0.000000 


« 

BOOT • 

SECTIONAL DERIVATIVES • 


MODE NO I 


NO. 

V 

1 

VOS 

lift COEFFICIENT 

MOMENT COEFFICIENT 






REAL 

IMAG 

REAL 

IMAG 

1 

2.0000 

-.5000 

2.0000 

0.000000 

0.000000 

0.000000 

0.000000 

2 

2.0000 

-.5000 

2.0000 

0.000000 

0.000000 

0.000000 

0.000000 

3 

0.0000 

0.0000 

0.0000 

0.000000 

0.000000 

0.000000 

0.000000 


TOTAL derivatives • 


C2 


FORCE COEFFICIENT IN I DIRECTION 
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CY • FOFCe COEFFICIEMT IN Y OIRECTION 
CM m PITCHING MOMENT COEFFICIEKT 
CN - YPwING MOMENT COEFFICIENT 
CL I • POLLING MOMENT COEFFICIENT 

PEAL I HAG 


Cl • 0.000000 0.000000 
CM • 0.000000 0.000000 
CL (POLL) « 0.000000 0.000000 


CENERALUEO AIRFORCES 


PRESSURE MODE DEFLECTION MODE 


1 

2 

3 


I 

1 

1 
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KOOE KO 2 


BOX PHESSURE DIFFERENCES 


BOX NO. 

XOC 

X 

V 

2 

PRESSURE difference 
REAL IMAGINARY 


1 

. 125C0 

2.35938 

1. IA950 

.50000 

.91A9A382E»01 

0. 

2 

• t 2500 

3.296R8 

1.1A950 

.50000 

.310775AlE»0l 

0. 

3 

.12500 

2.57813 

1. 71650 

.50000 

.113A36A0£»02 

0. 

A 

.62500 

3, 39063 

1.71650 

.50000 

.30A87371E*01 

0. 

S 

.12500 

2. 79688 

2.25000 

.62500 

.75669673E»01 

0. 

6 

.62500 

3.A8A36 

2.25000 

.62500 

. 1A97 1521E*0 1 

0. 

7 

.12500 

3.01563 

2.75000 

.87500 

.699'lil09lE*0l 

0. 

8 

.62500 

3.57813 

2.75000 

.67500 

.lO673323F»0l 

0. 

9 

.12500 

2.68750 

2. DOQOO 

.25000 

.5l312678E*0l 

0. 

10 

.62500 

3.A3750 

2.OOU0O 

.25000 

.15097670E*01 

0. 

1 1 

.12500 

2.68 750 

2.21650 

-.12500 

-.3326988AE*0l 

0. 

12 

.62500 

3.A3750 

2.21650 

-.12500 

. 1128A985E«01 

0. 

13 

. 1 2500 

2.68750 

2.A3300 

-.50000 

-.11A971A6E«00 

0. 

lA 

.62500 

3. A3 750 

2.A3300 

-.50000 

-.A0D36 1 1 7E»00 

0. 

IS 

. 12500 

2.68750 

2.21650 

-.87500 

.2A3l6808Et02 

0. 

16 

.62500 

3.A3750 

2.21650 

-.87500 

-.209690A3E»0 1 

0. 

17 

.12500 

2.68750 

1. 78350 

-.87500 

.15A0 7 71 1E4’01 

0. 

IS 

.62500 

3.A3750 

1. 78350 

-.87500 

-.23831 75AE*0 1 

0. 

19 

.12500 

2.68750 

1.56700 

-.50000 

-.3156A297E*01 

0. 

20 

.62500 

3.A3750 

1.56700 

-.50000 

-.M5071A9E»01 

0. 

21 

.12500 

2.68750 

1. 78350 

-.12500 

-.825761356*01 

0. 

22 

.6 2500 

3.A3750 

1. 78350 

-.12500 

-.27'327022E*00 

0. 

23 

.12500 

2.25000 

.21650 

.87500 

.201362626*01 

0. 

2A 

.62500 

3.25000 

.21650 

.87500 

.101623A6E*01 

0. 

25 

. 12500 

2.25000 

.6A950 

.62500 

.2A77881AC*01 

0. 

26 

.62500 

3.25000 

.6A950 

.62500 

.10735329E«01 

0. 

27 

.12500 

2.25000 

.86600 

.25000 

-.A2AA2539E*0l 

0. 

26 

.62500 

3.25000 

. 86600 

.25000 

-. 1956630AE*01 

0. 

29 

.12500 

2.25000 

.86600 

-.25000 

-.296AA657E *01 

0. 

30 

.62500 

3.25000 

.86600 

-.25000 

-.167838756*01 

0. 

31 

. 12500 

2.25000 

.6A950 

-.62500 

-.180667166*01 

0. 

32 

.62500 

3.25000 

.6A950 

-.62500 

-.89A02A65E*00 

0. 

33 

.12500 

2.25000 

.21650 

-.87500 

-. 1A8609A2E*01 

0. 

34 

.62500 

3.25000 

.21650 

-.87500 

-.669534666*00 

0. 



PANEL • 

SECTIONAL OEHIVATIVES • 


NODE NO 2 


STRIP NO. 

V 

1 

YOS 

1 

1.1A95 

.5000 

1.1A85 

2 

1.7165 

. 5000 

1.7165 

3 

2.2500 

.6250 

2.2500 

A 

2.7500 

. 8 750 

2.7500 

5 

2.0000 

.2500 

2.0000 

6 

2.2165 

-.1250 

2.2165 

1 

2.A330 

-.5000 

2.A330 

B 

2.2165 

-.6750 

2.2165 

9 

1p7B35 

-.6750 

1.7835 

10 

1.5670 

-.5000 

1.56 70 

11 

1. 7835 

-.1250 

1.7835 

12 

.2165 

. 8750 

.2165 

13 

.6A95 

.6250 

.6A95 

lA 

. 866 0 

.2500 

.8660 

15 

.6660 

-.2500 

.8660 

16 

.6A<J5 

-.6250 

.6A9S 

17 

.2165 

-.8750 

.2165 

• 

• 

BODY 

• 

• 


« sectional DERIVATIVES 

• 




# 






MODE 

NO 2 

BODY 

NO. 

Y 

2 

YOS 


1 

2.0000 

-.5000 

2.0000 


2 

2.0000 

-.5000 

2.0000 


3 

0.0000 

O.OOOO 

0.0000 

♦ 





• 

TOTAL 

DERIVATIVES 



• 



• 



C2 - FORCE COEFFICIENT IN / DIRECTION 


H* 


LIFT COEFFICIENT NOMENT COEFFICIENT 


REAL 

MAG 

PEAL 

I HAG 

6.128596 

0.000000 

-1.5A301 3 

0.000000 

7.196189 

o.nonooo 

-1.661 708 

0.000000 

A. 532060 

o.oououo 

-.9A0795 

0.000000 

A. 0335 7 1 

0.000000 

-.771029 

0.000000 

3.120518 

o.ocoono 

-.792507 

o.oooouo 

■1.0992A5 

0.000000 

-. 1AA719 

0.000000 

-.257666 

o.ocoooo 

.132299 

0.000000 

. 167368 

0-000000 

.503303 

0.000000 

-.A21202 

o.onoooo 

.6A8AAA 

0.000000 

-2.25 35 72 

o.ocoooo 

.619375 

0.000000 

■A.268AA2 

0.000000 

.603373 

0.000000 

1.5 1A930 

o.oocooo 

-.AA3A25 

o.oocooo 

1.775707 

0.000000 

-.A903A7 

0.000000 

•3. 100AA2 

0.000000 

.876713 

0.000000 

■2.321A27 

0.000000 

.709775 

0.000000 

■1.350 3A8 

0.000000 

.392300 

0.000000 

•1.077815 

0.000000 

.302111 

0.000000 


LIFT COEFFICIENT 
REAL MAG 

0.000000 0.000000 

.000000 0.000000 

0.000000 0.000000 


MOMENT COEFFICIENT 

PEAL I MAG 

0.00000 0 0.0 0 0000 

.163363 0.000000 

.11635$ 0.000000 




CV - force coefficient in ¥ DIRECTION 
CM ■ PITCHING moment coefficient 
CN • YAWING MOMENT COEFFICIENT 
CLI • ROLLING MOMENT COEFFICIENT 

real I mag 


Cl • 1502^2 0.000000 
CM • -A.2S)0<V0 0.000000 
CUFOLLi ■ 0.000000 0.000000 


GENERALIZED AIRFORCES 


PRESSURE MODE DEFLECTION MODE 


2 

2 

2 


1 

2 

3 


REAL IMAG 


CY • -.357263 0.000000 

CN - 0.000000 0.000000 


GENERALIZED FORCES 

REAL IMAGINARY 

-.201589A0E»02 0. 

-.61205394E*02 0. 

-.28271777E*02 0. 



I 


[ NO. 

xoc 

X 

Y 

1 

.12500 

2.35938 

1.14950 

2 

.62500 

3.29608 

1. 1S950 

3 

.12500 

2.57H13 

1.71650 

5 

.o2S00 

3.39063 

1. 71650 

5 

.12500 

2.79688 

2.25000 

6 

.62500 

3.4B43S 

2.25000 

7 

.12500 

3.01563 

2. 75000 

8 

.62500 

3.57813 

2. 75000 

9 

. 1 2500 

2.68750 

2.00000 

10 

.62500 

3.43750 

2.00000 

11 

.12500 

2.68750 

2.21650 

12 

.62500 

3.43750 

2.21650 

13 

.12500 

2.66750 

2.43300 

14 

.62500 

3.43750 

2.45300 

15 

. 12500 

2.68750 

2.21650 

16 

.62500 

3.43750 

2.21650 

17 

.12500 

2.61750 

1.78350 

IS 

.62500 

3.43750 

1.78350 

19 

.12500 

2.6-1750 

1.56700 

20 

.62500 

3.63750 

1.56700 

21 

.12500 

2.66750 

1.78350 

22 

.62500 

3. *3 750 

1.78350 

23 

.125CO 

2.25000 

.21650 

24 

.62500 

3.25000 

.21650 

25 

.12500 

2.25000 

.64950 

26 

.62500 

3.25000 

.64950 

27 

.12500 

2.25000 

.86600 

28 

.62500 

3.25000 

.86600 

29 

.12500 

2.25000 

.86600 

30 

.62500 

3.25000 

. 86600 

31 

*12500 

2.25000 

.64950 

32 

.6 25 00 

3.2500Q 

,64950 

33 

.12500 

2.25000 

.21650 

34 

.62500 

3.25000 

.21650 


BOX PRESSURE DIFFERENCES 


3 


NO 3 


PRESSURE DIFFERENCE 


REAL INAGINAAV 

.SOOOO 0. 0. 
.SOOOO 0. 0. 
.50000 0. 0. 
.50000 0. 0. 
.62500 0. 0, 
.62500 0. 0. 
.BTSOO 0. 0. 
.87500 0. 0. 
.25000 0. 0. 
.25000 0. 0. 
.12500 0. 0. 
'.12500 0. 0. 
.50000 0. 0. 
.50000 0. 0. 
.BTSOO D. 0. 

. BTSOO 0. 0. 
.87500 0. 0. 
.BTSOO 0. 0. 
.50000 0. 0. 
.50000 0. 0. 
.12500 0. 0. 
.12500 0. 0. 
.BTSOO 0. 0. 
.STSOO 0. 0. 
.62500 0. 0. 
.62500 0. 0. 
.25000 0. 0. 
.25000 0. 0. 
.25000 0. 0. 
.25000 0. 0. 
.62500 0. 0* 
.62500 Q. 0. 
-.87500 0. 0. 
-.87500 0. 0. 


h-* 

§ 


• 

• • 


• PANEL 

• 


« SECTIONAL DERIVATIVES 

• 


• 

• • 

MODE NO 3 

TRIP NO. Y 

1 

YOS 


1 

I.1AV5 

.5 000 

1 .1595 

2 

1.7165 

.5000 

1.7165 

3 

2.2500 

.6250 

2.2500 

A 

2.7500 

. 8 750 

2.7500 

s 

2.0300 

.2500 

2.0000 

t> 

2.2 165 

- . 1250 

2.2165 

7 

2.5330 

-.5000 

2.5 330 

8 

2.2165 

-.8750 

2.2165 

V 

1.7835 

-.8750 

1.78)5 

10 

1.5670 

-.5000 

1.5670 

1 1 

1.7835 

- . 1 250 

1.7835 

12 

.2165 

.8 750 

.2155 

1 3 

.65‘I5 

.6250 

.6595 

lA 

. 8660 

.2500 

.8660 

IS 

.8660 

-.2500 

.8660 

16 

.6595 

-.6250 

.6595 

17 

.2165 

-.8750 

.2165 


bOOY 

seciiokal derivatives 


HQOE KQ 3 


eODV NO. V 

2 

vOS 

1 2.0000 

-.5000 

2.0000 

2 2.0000 

-.5000 

2,0000 

1 0.0000 

0.0000 

0.0000 

• 

• total D£(i[VAnvES 

• 

• 

# 



Zl • FORCE COEFFICIENT In I OI«>ECIION 


LIFT Coefficient momcnt coefficient 


REAL 

[MAC 

heal 

1 H kC 

0.000000 

0.000000 

0. OOOOOO 

0.000000 

Q .OOOOOO 

0.000000 

0.000000 

0.000000 

O.OOOOOO 

0.000000 

0.000000 

0.000000 

0.000000 

o.oooooo 

o.ooooco 

0.000000 

0. OOOOOO 

0.000000 

o.oooooo 

0.000000 

o.oooooo 

o.oooooo 

0.000000 

0.000000 

O.OOOOOO 

o.oocooo 

0.000000 

0.000000 

o.oooooo 

o.oooooo 

0. OOOOOO 

o.oooooo 

0.000000 

o.oooooo 

0.000000 

o.oooooo 

0 .oOOOOO 

o.oooooo 

0.000000 

o.oooooo 

o.oooooo 

0 . OOOOOO 

o.oooooo 

0. OOOOOO 

0.000000 

0.000000 

0.000000 

o.oooooo 

o.oooooo 

0.000000 

o.oooooo 

0.000000 

0.000000 

0.000000 

0. OCOODO 

0.000000 

o.oooooo 

0.000000 

0.000000 

o.oooooo 

o.oooooo 

o.oooooo 

0.000000 

0.000000 

O.oooooo 

O.oooooo 

0.000000 

0.000000 


LIFT COEFFICIENT MOMENT COEFFICIfNT 

PEAL IMAC Ft*C 


0.000000 

0.000000 

0.000000 


0.000000 

o.oocooo 

O.ODOODO 


0.000000 
0. 000000 
0.000000 


0.0 00000 
0.000000 
0.000000 



i 

i 

i 

1 

I 

I . , . 

( 

! CY « FOI^CC CUEFFtClCHT IH Y 01 DECT 1 OH 

1 CM - PirCHiNC HOHEMT COEFFICIENT 

CN > YAHINC MONENr COEFFICIENT 

j cii ■ polling moment coefficient 

' REAL. IN«C 


C2 « 0. 000000 0.000000 

CM m 0.000000 0.000000 

CUFOLLI - 0.000000 0.000000 


• GCIIERALItED AinFORCES * 


PRESSURE node deflection M 0 D£ 



REAL (HAG 


CY - 0.000000 0.000000 

CN » 0.000000 0.000000 


CENERAtlZEO FORCES 
REAL IhaCINARV 



r 


o 


* 

• PEOUCCD FPCOUEHCV ■ .fO000(.*00 


0. 

THE )40DUNUASH 
. 66666 TE *00 0. 

W ELEHENTS for mode no. 1 

.6666678*00 0. 

,666667E*00 

0. 

.666667E *00 

0. 

,S96000E*00 

0. 

.S96000E*00 0. 

.S96000E»00 

0. 

.S96000E*00 

0. 

0. 

0. 

0. 

0. 

0. 0. 

0. 

0. 

O. 

0. 

0. 

0. 

0. 

0. 

0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

0. 

0. 0. 

0. 

0. 

0. 

0. 

0. 




THE 340at.NUASH U elements FQp hOOE no. 2 

.100000E*Ol .lBSS42E4ai .100OOOE*Ol .2%10«2E*01 .lOOOOOEtOT .19B<)Se£»01 .10000oE*01 .2S)12SE«01 .S99000E»00 .1$T1B1E»01 


.S94000£»00 .22SISbE«Ol .e94000£*00 .196494E»01 .S94000E«00 .210019E401 0. O. 0. 0. 


0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 




THE }4D0UNt<6SH 

U 

elements FCe HOOE NO. 

5 







0. 

.76633)E*00 0. 


. 766}3)E*00 0. 



. tl94J}E*01 

0, 

.114L33E*01 

0. 

.1S2767E*01 

0. 

.IS2767E*01 0. 


.1900}SE*01 0. 



.190035£*01 

0. 

1666678*00 

0. 

-.1666678*00 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0, 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 




IME 

9 

Hj iutpix Elements fop 

MODE 

NO. 

1 





0. 

0. 0. 


0. 0. 



0. 

0. 

.6C666TE*00 

0. 

.6666678*00 

0. 

. 66666 7E *00 0. 


•666667E*00 0. 



.6666&7E*00 

0. 

.6666678*00 




THE 

9 

MJ-PPIME elements fob 

MODE 

NO. 

1 





0. 

0. 0. 


0. 0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 0. 


0. 0. 



0. 

0. 

0. 



COL. NO. 

1 OF LINE ELEMENT FPE5SUPE COEFFICIENT HtTFIX 







0. 

0. 0. 


0. 0. 



0. 

-.349066£*00 

.4188798*01 

-.6991328*00 

0. 


•,}4«0t«E»60 '«16?5S2£»91 -.69S132C*00 *2099V0£»(H 1 J9l.2frf «01 0. -.<>991 J2E*00 -*2O9940ErOl 


THE 9 HJ HATKIX elements FOP MODE NO. 


2 



0. 


0. 0* 

0* 0. 


0. 

.100000F«01 

*1$OOOOE«01 

.100000F*01 

.19166T£*01 

.IOOOOOE401 

.2S8333E40L .lOOOOOE^OL .bB8b&TC400 »lOOOOaE«aL 

.Za000QE*01 

.lOOOOOEvOl 

♦333333E *01 





THE 9 

hj-ppihe elements rop mode 

NO. 

2 





0. 


0. 0. 

0. 0. 


0* 

0. 

.6666675*00 

0. 

.666667E*00 

0. 


.t66t67E»00 0. 

,666667E»00 0. 


.6666&T£»00 

0. 

.6666675*00 



COL. hO. 

2 OF LINE ElEhENT PFESSuHE COEFFICIENT HATMX 







0. 


0. 0. 

0. 0. 


0. 

.5«97T9E*01 

.109720E*02 

-.200U3E*01 

.209660E*01 

-.3B6590E»OI 

-.54A)b3E»0L .24b3C6E401 .AlSfiT9E*01 '.blSB79E*01 

.«iee79E»oi 

-*6632255*01 

-,837758E*01 





THE 9 

HJ HATBiX ELEMENTS FOB MODE 

NO. 

3 





0. 


.333333E400 0. 

*333333E»0D 0. 


.333333E*00 

0. 

.1333335*01 

0. 

.133333E*01 

0. 


.)3J3J3E*01 0. 

0. 0* 


0 . 

0. 

0. 





THE 9 

WJ^PftlME ELEMENTS FOR MODE 

HOm 

3 





0. 


0. 0. 

0. 0. 


0* 

0. 

0. 

0. 

0* 

0. 


0* 0. 

0. 0. 


0. 

0. 

0. 



CQL« HQ* 

3 OF LINE ELEHENT PltE»U«E COEFFICIENT HfcTftlK 







-^.174S33£400 

.209A40E401 ~.34906A£400 0. «.IT4S33E*00 

-.83775BE400 

-.6901325*00 

.a37rsoE*oi 

-.1396265*01 

0. 

-.b9aX32E*00 

-.33S103E401 0. 

0. 0. 


0. 

0. 

0. 




PANtL 

SECriOUAL 0€A[tr4TIVES 


MODE HO 1 


SlRIA NO. 

V 

2 

Y05 

Llfl COEfFlClENl 

KOHfNT COEFFICIENT 





REAL 

i^tr. 

«FAl 

MAC 

1 

1. lAOJ 

.5000 

1.15Y5 

-1 .0 7 7532 

3. l«490 1 

. 790016 

-1.015790 

2 

1. M(>$ 

.5000 

1.7165 

-1 . 32 701 7 

3.704252 

.891 397 

-1 .001647 

3 

7.2iOO 

.6250 

2.2500 

-.652704 

2.570390 

.500485 

-.710637 


2. 7BC0 

.8 750 

2.7500 

-.280565 

7.235104 

.353334 

-.561 7J2 

5 

2.0000 

.2 500 

2.0000 

-.52H1 7 

1 .4)11616 

.303904 

-.410516 

6 

2.71GB 

-.1250 

2.2165 

). 199165 

- .567 326 

-.739911 

-. J5"466 

7 

7. ‘.330 

-.5000 

2.6330 

.9M367 

-.143000 

-.201942 

-. lOOTM 

B 

2.2165 

-.8750 

2.2165 

-1 . 1059 70 

. 107618 

- 1 

.175621 

•} 

1.7«39 

-.0750 

1.7015 

-1.101796 

'.C60633 

.oontif' 

.226432 

10 

1.5670 

-.5000 

1.5470 

1.122871 

• .929562 

-.449241 

. 1 14996 

1 1 

1.7035 

-.1250 

1.7P 55 

3.600597 

-l.*; 72065 

- 1 .104455 

.0140..9 

12 

.2165 

.8 750 

.2165 

.640020 

-.00451 7 

- , 341064 

-.132291 

13 

.65<»5 

.6250 

.6665 

.520005 

.212203 

- . 264279 

-.186194 

I* 

.0660 

.2 500 


[.122890 

-2.1 16290 

7506*4 

.650997 

1) 

.6660 

-.2500 

.6660 

.651506 

-1.56 7688 


• 52 7301 

!« 

.6*95 

-.6250 

.6405 

.095301 

- .020265 

-.149001 

.293226 

IT 

.2165 

-.8750 

.2165 

.005656 

-.573246 

-.092716 

.199)00 

• 

• 

BODY 

• 






• SECTIOHAL DEKlVAtlveS 

♦ 

* 

• 









NODE 

NO 1 





BODY NO. 

Y 

Z 

V05 

lift coefficient 

HQHENI COEFFIClENt 






real 

IKAG 

«EAL 

t NAG 

1 

2.0000 

-.5000 

2.0000 

0.000000 

0.000000 

0.000000 0. 

000000 

2 

2.0000 

-.5000 

2-0000 

-.132645 

-.000000 

.056090 

100909 

3 

0.0000 

0. 0000 

0.0000 

-.237711 

0*000000 

.116355 

0775 70 


lOiAi OCeivjirjVEJ 


Cl - fOACE COEFFJCJfWT jfl I DIAfCTIDN 



I 


I cv • FOBCC COEFFICICNT in Y OIPECTION 

1 CN . PIlCHlNC NOHENI COEFFICIEM 

1 CN - VANtNC HOHENT CDEFFICIEM 

I Cl I > ROILINC HOHENr COfFFIClENf 

I PEAL IHAC 


C/ • -.16709* 1.991560 
CH - .45*503 -1.967*75 
CL (ROLL) • 0.000000 0.000000 


« GtNERALUEO AIRFORCES • 



I 



PRESSURE MODE DEFLECTION NODE 



REAL (HAG 


CV .312622 -.5*6109 

CN > D.OODODO 0.000000 


GENERALIZED FORCES | 

REAL IHACINARV 

.107*3960E*01 -.95**67115*01 1 

.62806L27E*01 - . L0fla33l5E *02 t 

-.5207302*E*01 15261 *35E*02 I 

i 


I 


114 


PAN^L • 

JfCMC’NAl. ofpiwnves • 


HQOl na 2 


SIMP NO. 

Y 

7 

YOS 

1 

1.1*95 

.5000 

1. I *95 

2 

1. U*5 

. 5000 

1.7165 

1 

7 . 7 500 

.6250 

2.2500 

4 

7. 7500 

.0 750 

2. 7520 

5 

2.0000 

-2 500 

2,0000 

6 

2.2 165 

-.1250 

2.2165 

7 

7.5130 

-.5000 

2. <.310 

8 

2.2165 

-.8150 

7.2165 

9 

1. 7flJS 

-.0750 

1.7015 

10 

1.5670 

-.5000 

1.5670 

ll 

1.7635 

-. 1250 

1.7835 

17 

.7165 

. 8750 

.2165 

13 

.6595 

.6250 

.6*95 

1* 

. 8660 

.2 500 

.8660 

15 

. 8660 

-.2500 

.8660 

1* 

.6*95 

-.6250 

.6*95 

17 .71*5 

♦ 

♦ Root 

< s ec r fOMA L oc^ r VA r c ve s 

• 

-.8>50 

• 

• 

.21*5 




HOOl 

NO 2 

aooY HO- 

t 

2 

ros 

1 

2.0000 

-.5000 

2.0000 

2 

2.0000 

-.5000 

2. 0000 

3 

o.oooo 

0.0000 

0.0000 

* 


♦ 


* total 

0€»I VATIVES 



♦ 


• 



cz - fOt»Cl COlfriCl^HZ In Z OfPfcUON 


Lin coEFficitNi HOHfii cocfucifni 


BEAL 

)-*C 

BF*L 

1 kOC 

1 .990107 

12. 57*1 77 

. 5 7666* 

-5. 155086 

1 .9 JI919 

1*. 591,729 

. 700990 

-5.522tjt, 

2.05*5*5 

9.BM698 

. 3 l*f 09 

-3.532105 

2.5* n JO 

B.222CDI 

.1 92291 

-2, 752 M 2 

.059373 

5. M0992 

.tn* 7 75 

-2.12*639 

9 .* 7 IT I ; 

- 7 . *92 U 1 

-3.1 12902 

-.1 71.690 

2 .0*1166 

-2 , *0006* 

-1 .16OC07 

..’9" 221 

-3.63*591 

1 .053*69 

. 6*?567 

.962125 

-3.91232? 

1 .29 7.1*0 

.*7*15* 

I . 1 709M 

2.090561 

-5.2010*3 

-1.295902 

1.3551*7 

0.613097 

-12.91*563 

-3.552102 

1. 79 7665 

1 .90*670 

-1.1*6360 

- 1 . >59560 

.167016 

1.9929** 

-.215539 

-5.227T1 1 

-.20**09 

-.006960 

-a. 6**97* 

-1,200906 

3.5776*1 

-.566900 

-6.102*62 

-.631969 

2.627053 

-1.061793 

-2,71*130 

.00937* 

1.219317 

-.750193 

-5-7*7225 

,070337 

. 757299 


LIFT COfFFJCIfKI COrFFICIFNT 

PEAL IHAC “F*l 


0.000000 

-, A 0 Tt 3 A 

-.49SJ3? 


O.OOOODO 

. W ' l/AA 

I 


0.000000 
. M 7135 


0.000000 
.: f B555 
. I 1*355 



CV - FOfcCE COEFFICIENT IN Y OIPECTIDN 
CH - PITCHING HOMENl CnEFFICIENT 
CN • YAWING MOMENT COEFFICIENT 
CLI > FGLLING moment COEFFICIENT 



REAL 

I NAG 


REAL 

IHAG 

CZ - 

Z. 0^915A 

A. TA68AS 

CV - 

.099798 

>2.336862 

CM • 

-UZ20O99 

-7.33S9B2 

CN ■ 

0.000000 

O.OOOOQQ 

CLIRQLLI 

• 0.000000 

0.000000 





GENEHALIZEO AIRFORCES 


PRESSURE MODE 

DEFLECTION MODE 

GENERALIZED 

FORCES 



REAL 

IMAGINARY 

2 

1 

-. I3A32TAAEP02 

-.303TA63AE»02 

2 

2 

>.U70SOOOE»02 

-.70A0a7TlE*02 

2 

3 

-.422190A9EA02 

-.4244612ZEA02 
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